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THE SPECIFIC HEATS OF CRYSTALS 
Part I. General Theory 
By Sirk C. V. RAMAN 


(Memoir No. 87 of the Raman Research Institute, Bangalore 6) 
Received October 12, 1956 


1. INTRODUCTION 


THE quantum theory of the specific heats of crystals proposed by Einstein 
in the year 1907 regards a crystal as an assembly of an immense number of 
harmonic oscillators with specific frequencies, the energy of vibration of 
which obeys the quantum rule; these oscillators form a system in thermo- 
dynamic equilibrium of which the behaviour can be described statistically 
with the aid of Boltzmann’s theorem. Einstein identified these oscillators 
with the structural units (“‘ Elementargebilde ”’) in the crystal and assumed 
that they could be grouped into sets, each set comprising a great number of 
oscillators characterised by a common frequency of vibration; the total 
number of oscillators of all sorts is taken to be thrice the number of atoms 
comprised in the crystal. He showed it to be a consequence of his theory 


that the thermal properties of a crystal would stand in the closest relationship 
with its spectroscopic behaviour. 


The assumptions on which Einstein’s theory rests are clearly justified 
by the physical facts of the case. For, every crystal consists of an immense 
number of structural units similar to each other which are capable of mecha- 
nical vibration and which by reason of their similarity may be expected to 
exhibit identical dynamical behaviour. The theory is, in effect, a synthesis 
of the results indicated by classical dynamics for such a system with the 
notions of the quantum theory and the basic principles of thermodynamics. 
It involves no inherent contradictions and can therefore claim to be a rational 
approach to the fundamental problems of crystal physics. The following 
questions, however, need to be answered before any use can be made of the 
theory. How are the oscillators of different sorts envisaged in the theory 
telated to the known structure of a crystal ? How are they to be enumerated ? 
What is the procedure by which the characteristic frequency of each set of 
oscillators may be evaluated ? The present memoir concerns itself with finding 
the answers to these questions and thus giving form and substance to the 


ideas of a very general nature put forward by Einstein. 
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2. THE NORMAL VIBRATIONS OF CRYSTAL STRUCTURES 


That the specific heat of a crystalline solid tends with rising temperature 
to reach a limiting value proportional to the number of atoms comprised in 
it is a clear indication that the theoretical approach to any explanation of the 
facts has of necessity to consider the subject from an atomistic standpoint, 
We have also to recognize the existence of interatomic forces; hence, the 
movements with which we are concerned are not simple translations of the 
individual atoms but vibrations of the atoms about their positions of equi- 
librium in the structure of the crystal. Any attempt to describe these vibra- 
tions or to enumerate them has necessarily to be based on the fundamental 
theorem in classical mechanics which states that all the possible small vibra- 
tions of a connected system of particles about their positions of equilibrium 
are summations of an enumerable set of normal modes; in each such normal 
mode, the particles of the system have all the same frequency and the same 
or opposite phases of vibration. Hence, any rational approach to the specific 
heat problem must of necessity start with a consideration of the normal 
modes of vibration—in the sense of the theorem just stated—of the atoms 
about their positions of equilibrium in the lattice structure of the crystal. 


The basic principle of crystal architecture is that the structure comes 
into coincidence with itself following a unit translation along any one of the 
three axes of the lattice. It follows as a necessary consequence that the 
normal modes of vibration of the atoms characteristic of the structure 
of the crystal should possess the same property. This can evidently 
happen in two ways, thereby enabling us to divide the normal modes 
into two distinct classes. In the first class of normal modes, the 
amplitudes as well as the phases of oscillation of the equivalent atoms 
which come into coincidence following a unit translation are identical. 
In the second class of normal modes, the amplitudes of equivalent 
atoms are the same but the phases are all reversed and the normal 
mode therefore remains the same following the unit translation. Since these 
two alternatives are possible for a unit translation along each of the three 
axes of the lattice, we have 2 x 2 x 2 or 8 possible situations. Consider- 
ing each of these situations separately, we proceed to write down the equations 
of motion of the p atoms comprised in the unit cell in terms of the 3 p co- 
ordinates which determine their displacements. Since the displacements of 
the equivalent atoms in adjoining cells which interact with the atoms in the 
cell under consideration are also the same, the 3 p equations of motion of 
the latter contain only 3p independent co-ordinates. They can therefore 
be completely solved, giving us 3 p solutions for each of the eight possible 
situations referred to above. Thus in all, we have 24 p solutions. (3 p —3) 
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of these solutions represent normal modes of vibration of the first kind, 21 p - 
solutions represent normal modes of the second kind and the three remaining 
solutions represent simple translations which have necessarily to be 
excluded. 


The (24 p — 3) normal modes of vibration of the atoms in a crystal 
indicated by the theory may be regarded as the modes of internal vibration 
of the group of 8 p atoms comprised in a super-cell of the crystal lattice whose 
linear dimensions are twice as large as that of the unit cell consisting of 
patoms. The three omitted solutions would then represent the translations 
along each of the three axes of the crystal lattice of the whole group of 8 p 
atoms included in the super-cell. Thus, we recognise that the structural 
unit whose dynamical behaviour is representative of the entire crystal is not 
the unit cell of the crystal structure, but has twice its dimensions along each 
of the three axes of the lattice. 


Thus, following Einstein, we may regard a crystal as an assembly of 
(24 p — 3) sets of oscillators, each of the sets having its own distinctive fre- 
quency of vibration. The number of oscillators in each set is the same as 
the number of super-cells containing 8 p atoms which are included in the 
whole crystal. Each set of oscillators forms an assembly in thermodynamic 
equilibrium and its behaviour can therefore be described statistically with 
the aid of Boltzmann’s theorem. The average energy of an oscillator in each 
set is given by Einstein’s formula for the particular frequency. Multiplying 
this by the number of super-cells contained in the volume of the crystal, we 
obtain an expression for the energy of all the oscillators in that set. Finally, 
summing up the expressions thus obtained for all the (24 p — 3) sets with 
their respective frequencies of vibration, we obtain an expression for the 
thermal energy arising from the excitation of the vibrations of the atoms having 
precisely specifiable frequencies. 


3. THE VIBRATION SPECTRA OF CRYSTALS 


The foregoing results may be summed up as follows. If a unit volume 
of the crystal contains S super-cells each including 8 p atoms, the total number 
of degrees of atomic freedom is 24 pS; a very large proportion of this number, 
viz., (24 p — 3)S degrees of freedom, appears in the vibration spectrum of 
the crystal as sharply defined monochromatic frequencies; (3 p — 3) of these 
frequencies represent vibrations in which equivalent atoms contained in each 
super-cell have the same amplitudes and phases of vibration, while the 21 p. 
other frequencies represent modes of vibration in which the amplitudes are 
the same but the phases are opposite for equivalent atoms along one, two or 
all three of the axes of the lattice. 
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We have now to consider and determine the nature of the movements 
represented by the 3S degrees of atomic freedom left out in the preceding 
enumeration. This is also the number of degrees of freedom of translation 
of the S super-cells included in the crystal. Any translation of an individual 
super-cell would necessarily set up forces resisting such movement and also 
tending to displace neighbouring super-cells. We are therefore led to assume 
that the 3S degrees of freedom under discussion would manifest themselves 
as internal modes of vibration in elements of volume of larger dimensions 
than our super-cell. The larger such element of volume is, the lower would 
be the limit of the possible frequencies of such vibration. The translations 
of the individual super-cells are therefore to be associated with a whole series 
of frequencies lying between zero and an upper limit which is set by the fre- 
quency of internal vibrations in our super-cell. To take account of the ani- 
sotropy of the crystal, we would, in general, have to assume separate upper 
limits of frequency for the translations along the three axes of the crystal, 
It is left to determine how the 3 S unassigned degrees of freedom are distri- 
buted amongst the oscillators of various frequencies lying in the permitted 
ranges O-—>vy, O— vy, O— vz. 


The problem indicated above can be dealt with and solved in a very 
simple manner. Consider two volume elements in the crystal which have the 
same shape as our super-cell and whose dimensions are respectively m and n 
times larger. The number of degrees of freedom of translation of these en- 
larged super-cells is respectively three times the number of such cells included 
in the volume of the crystal; the difference between them, viz., 3 S(1/m' 
— 1/n*) would therefore be the number of oscillators whose frequencies lie 
between vm and vy, which are respectively their lowest frequencies of inter- 
nal vibration. But the principle of dynamical similarity enables us to write 
vm = k/m and vy = k/n where k is a constant having the dimensions of a 
frequency. Hence, as m approaches n, vm — vpn may be set equal to dv and 
the number of oscillators whose frequencies lie between vm and vy, may be 
written as 9 S v? dv/k*. On integrating this between zero and v; (the upper 
limit of the permitted frequencies), we should recover 3 S as the total number 
of oscillators, and hence k® is identified as v;3. Accordingly, the law of dis- 
tribution of frequency in the spectrum may be written as 9S. v2dv/v,° for 
the isotropic case. More generally, for an anisotropic crystal vj; may be 
replaced by vy or vy or vz respectively for the three directions of translations 
and the numerical factor 9 is replaced by 3. 


The law of distribution of the frequencies thus derived indicates that 
the majority of the oscillators now under consideration have frequencies not 
far removed from the upper limit. The same result may be expressed other- 





















ce, a 


oOo 


tl 














The Specific Heats of Crystals—I 157 


wise by the statement that the majority of the oscillators are of very small 
dimensions not far removed from the lower limit, namely the dimensions 
of the super-cell. If, for example, we were to identify a// the vibrations arising 
from the translations of the super-cell with the internal vibrations in a volume 
element whose dimensions are twice those of the super-cell, we would not 
be seriously in error. All except one-eighth of the 3S degrees of freedom 
would be accounted for in such an enumeration; the residue omitted would 
represent oscillations in still larger volume elements appearing with still lower 
frequencies. 


4. THe THERMAL ENERGY OF CRYSTALS 


The atomistic approach to the specific heat problem proposed by Einstein, 
when fully developed, thus leads us to the following picture of the nature 
of the vibrations in a crystal which are the carriers of the thermal energy. 
The oscillating units are not the unit cells of the structure, but are twice as 
large as the unit cells in their linear dimensions; their oscillations are of two 
kinds, namely, those arising from their internal vibrations and those arising 
from their translatory movements. The oscillations of the former description 
appear with (24 p — 3) discrete frequencies of vibration, p being the number 
of atoms contained in the unit cells of the crystal structure. The translatory 
movements, on the other hand, appear with a spectrum of frequencies all of 
which are lower than those of the internal vibrations. These, for the most 
part, are concentrated near the upper limit of their permitted range; but 
there is a residue which tails off as a continuous spectrum to very low fre- 
quencies. 


Accordingly, making use of the reasoning employed by Einstein for 
evaluating the average energy of an oscillator in each of the sets under con- 
sideration, we obtain the following expression for the heat content of a volume 
of a crystal containing N unit cells of the crystal structure as a function of 
the temperature, namely 


ni DO 3 3hv*dv 

g ewakT —{ yp ) eeatT—y{ 
i=1 0 

The numerical factor 1/8 appears in the expression because each dynamic 


unit contains 8 unit cells of the crystal structure. On differentiating the 
expression with respect to T as usual, we get the formula for the specific heat. 


The following remarks may be made regarding the numerical evaluation 
of the expression given above. All the (24 p — 3) frequencies would be dis- 
tinct from each other only in the case of a completely anisotropic crystal, 
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If any symmetry elements are present, the number of distinct frequencies 
would be naturally diminished and the formula would then contain only such 
smaller number, but the individual terms would then have to be multiplied 
by their respective degeneracies. Such reduction in the number of distinct 
frequencies would naturally appear both in respect of the (3 p — 3) frequencies 
in which both the amplitudes and the phases are the same in adjacent unit 
cells and in the 21 p others in which the phases may be opposite as already 
explained. As an illustration of these remarks, we may consider a case in 
which p = 2 and the crystal belongs to the cubic class and its structure con- 
sists of two interpenetrating face-centred cubic lattices. The (3 p — 3) dis- 
tinct frequencies then reduce to a single triply degenerate frequency; the 
21 p or 42 other vibrations reduce to only eight distinct frequencies, all of 
which may be readily described in geometric terms related to the structure 
of the crystal. 


Since the first term in the expression for the thermal energy is a summa- 
tion extended over (24 p — 3) distinct modes of vibration, while the second 
represents only the three residual translations, it is evident that the latter 
would be of very minor importance relatively to the first, especially in those 
cases where p is very large, in other words when the crystal has a multi-atomic 
structure. The position is a little different when p is small, as for example, 
when p is equal to 1. The contribution from the second term would not then 
be altogether negligible in comparison with the first term. The second term 
also acquires some importance relatively to the first at very low temperatures. 
For, since the frequencies appearing in it are low, their contributions to the 
specific heat would survive when those due to the vibrations of the higher 
frequencies appearing in the first term have been effectively frozen out by 
reason of Boltzmann’s theorem. 


5. THE SPECTROSCOPIC BEHAVIOUR OF CRYSTALS 


In his paper of 1907, Einstein emphasised the intimate relationship 
indicated by his theory between the thermal properties of crystal and their 
optical behaviour. It is appropriate therefore that we make a reference 
here to the decisive support given to the results of the theory set forth above 
by the experimental facts regarding the spectroscopic behaviour of crystals 
revealed by investigations made with a great variety of materials and by 
diverse techniques of investigation. We need not however dilate on this topic, 
since the experimental situation in this respect has been reviewed in my 
Lindau address, and since, moreover, the particular case of the diamond which 
illustrates the general principles in a very striking fashion has been dealt 
with in detail in a recent lecture published in these Proceedings. 
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6. SUMMARY 


In Einstein’s atomistic approach to the specific heat problem, a crystal 
is regarded as an assembly of sets of individual harmonic oscillators, each 
set consisting of a great number of oscillators having a common frequency 
of vibration. In the present paper, this view is developed and an expression 
is obtained for the thermal energy of the crystal which appears as a summation 
of (24 p — 3) terms arising from the individual monochromatic frequencies 
of internal vibration of the dynamic units of the crystal structure each con- 
taining 8 p atoms, with a supplementary term which takes account of the 
translational movements of the same units manifesting themselves as a con- 
tinuous spectrum of frequencies of vibration within the crystal. 
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The Case of Diamond 
By Si C. V. RAMAN 


(Memoir No. 87 of the Raman Research Institute, Bangalore 6) 
Received October 15, 1956 


IN a lecture on the diamond published in the September issue of these 
Proceedings (Vol. 44, 1956, p. 99), the vibrations of the diamond structure 
were discussed in detail and the results of spectroscopic investigations on 
the nature of the vibration spectrum were also set forth. The results con- 
tained in that paper furnish us with all the material necessary for a complete 
theoretical evaluation of the specific heat of diamond over the entire range 
of temperatures for which experimental data are available. We shall now 
proceed to make such an evaluation and compare it with the results of the 
published measurements on the specific heat of diamond. 


Table I gives a list of the nine characteristic frequencies of the diamond 
structure expressed in wave-numbers (cm.~!), their respective degeneracies 
and the geometric description of the corresponding modes of vibration. The 
frequencies shown in the last column of the table were evaluated with the aid 
of the equations of motion and a set of force constants and were also inde- 
pendently confirmed by the results of spectroscopic investigation by several 


TABLE I 
Vibration Spectrum of Diamond 











No Degeneraey Description of the mode F oe in 

Es | 3 Oscillation of the two interpenetrating lattices against each | 1332 
other 

II | 8 Tangential oscillation of the octahedral planes | 1250 
Ill 6 Tangential oscillation of the cubic planes 1239 
IV | Normal oscillation of the octahedral planes 1149 
V&VI | (3+3) Normal oscillation of the cubic planes 1088 
Vil | 7 Normal oscillation of the octahedra] planes 1008 
Vill | 6 Tangential oscillation of the cubic planes 740 
IX | 8 Tangential oscillation of the octahedral planes 621 
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different techniques. It is a simple matter with the aid of the table to cal- 
culate the contributions to the specific heat arising from the nine modes of 
internal vibration of the diamond structure. The lowest frequency appearing 
in the table expressed in wave-numbers, viz., 621 cm.-! may also be properly 
assumed to be the upper limit of frequency in the continuous spectrum of 
the vibrations arising from the translations of the super-cell of the crystal 
structure, as explained in Part I of the present paper. 


The contribution to the specific heat arising from the continuous part 
of the vibration spectrum can be readily evaluated on this basis with the aid 
of the tabulated values of the relevant functions available in “ Landolt- 
Bornstein”. Table II below shows the results of such computations. The 
contributions to the atomic heat arising respectively from the discrete spec- 
trum and from the continuous spectrum are shown separately for a series of 
temperatures ranging from 15° to 1100° on the absolute scale. The sum of 
the two contributions is also tabulated. 


TABLE II 


Theoretical Computation of the Atomic Heat of Diamond 
(Calories per gram-atom per degree) 























| | | 
|Absolate | | Absolute | | 
| Tempe- | Discrete Continuous Tempe-| Discrete | Continuous 
| rature | Frequencies Spectrum Total rature |Frequencies | Spectrum Total 
Tv _ | \ 
| | 
15 ye 0-00014 0-00014 206 0-4069 | 0-1631 0-5700 
20 ‘ 0-00033 | 0-00083 | 225 | 0-5894 | 0+1895 | 0.7789 
25 me 0-00064 0-00064 250 0-7844 | 0-2124 0-9968 
30 os 0-0011 0-0011 275 1-0137 | 0-2319 1-2456 
40 os 0-0026 0-0026 300 1-2418 | 0-2486 1-4904 
50 .* 0-0051 | 0-0051 350 1-700 | 0-275 1-975 } 
60 a 0-0088 | 0-0088 400 2-139 0-294 2-433 
70 0-0005 0-0140 | 0-0145 450 2-538 | 0-308 2-846 
80 0-002 0-0207 | 0-0227 500 2-892 0-319 3-211 
90 0-0056 0-0290 | 00-0346 | 600 3-469 0-334 3+803 
| 100 0-0128 0-0388 | 00-0516 | 700 3-899 0-344 4-243 
| 120 0-0419 0-0619 | 0-1038 | 800 4-224 | 0-350 4-574 
| 140 0-0957 0-0877 0-1834 | 900 4-465 | 0-355 4-820 | 
| 160 0-1758 0-114 | 0-2898 | 1000 4-653 | 0-358 5-011 | 
| 175 0+2517 0-1333 | 0-3850 | 1100 4-797 | 0-360 5-157 | 
| | 











The values of the atomic heats of diamond listed in Table II are exhi- 
bited in graphical form in Fig. 1. At temperatures below 70° absolute, the 
curve appears to touch the axis since the very low values of the specific heats 
in this region cannot be represented on the scale of the graph. The values 
in this region are, however, of great interest, since, as will be seen from 
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Fic. 1. The atomic heats of diamond calculated from the spectroscopic data. 


Table I], they arise exclusively from the contribution of the continuous 
spectrum of frequencies which represents only 3 out of the 48 degrees of atomic 
freedom of movement of the crystal. Accordingly, the theoretical curve 
for this part of the temperature range has been drawn with a greatly enlarged 
scale for the ordinates in Fig. 2 below. Fortunately also, the specific heat 
of diamond in this region of temperature has been the subject of detailed 
investigation by Warren DeSorbo (Journal of Chemical Physics, Vol. 21, 
1953, p. 876) and his experimental results for this very interesting region are 
plotted in the diagram. A gratifying agreement between theory and observa- 
tion will be noticed. 


The specific heat of diamond in the range of absolute temperatures from 
70° to 300° was determined by K. S. Pitzer (Journal of Chemical Physics, 
Vol. 6, 1938, p. 68). The measurements by DeSorbo also cover this range. 
A careful inter-comparison between them shows that Pitzer’s values are dis- 
tinctly higher than those of DeSorbo throughout this range and that they 
also exhibit a scatter which is distinctly larger than the results of DeSorbo. 
DeSorbo’s work is of later date and his data form a continuous sequence 
over a wider range of temperatures. For these reasons, it was felt that a 
comparison between theory and experiment would be more appropriately 
made with the results of DeSorbo. Fig. 3 below shows such comparison and 
here again a very satisfactory agreement emerges. 










The Specific Heats of Crystals—II 





— THEORETICAL CURVE 
© OBSERVATIONS BY DESORBO 


20F 


/ DEGREE) 


3 
° 
7 


60F 


40k 


SPECIFIC HEAT X10 (GAL/Q7m. ATOM 


20r 














i i i 
20 40 60 
ABSOLUTE TEMPERATURE 





Fic. 2. Comparison between theory and experiment at very low temperatures. 
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Fic. 3. Comparison of theory and experiment below room temperature. 














164 Sir C. V. RAMAN 


Finally, Fig. 4 exhibits a comparison between the theoretical values and 
the experimental data in the temperature range from 300° to 1100° absolute 
covered by the determinations of Magnus and Hodler (Annalen Der Physik, 
Vol. 80, 1926, p. 808). The observed discrepancies throughout the range 
do not exceed two per cent. and in most cases are much less. The agreement 
should therefore be considered as satisfactory. 
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Fic. 4. Comparison between theory and experiment at high temperatures. 





SUMMARY 


The specific heat of diamond in the temperature range between 15° 
and 1100° absolute has been computed theoretically, making use of the nine 
spectroscopically determined vibration frequencies of the diamond structure 
as well as their known degeneracies, and further taking the lowest of these 
nine frequencies to be the upper limit of frequency in the residual con- 
tinuous spectrum as explained in Part I of the paper. The results are com- 
pared with the experimental observations of DeSorbo from 15° to 300° ab- 
solute and of Magnus and Hodler from 300° to 1100° absolute. A highly 
satisfactory agreement between theory and experiment emerges. 

















ON THE DIFFRACTION OF LIGHT BY 
PROGRESSIVE SUPERSONIC WAVES 


Oblique Incidence: Intensities in the Neighbourhood of the Bragg Angle 
By P. PHARISEAU§ 
(Rijksuniversiteit Gent, Seminarie voor Wiskundige Natuurkunde, Ghent, Belgium) 


(Communicated by Sir C. V. Raman, F.R.S., N.L.) 


1. INTRODUCTION 


IN the case of oblique incidence, Bhatia and Noble! have shown that, if the 
angle of incidence ¢% is close to the Bragg angle for the order —1, the 
expressions for the intensities differ from those in the case that the angle of 
incidence is very different from the Bragg angle. These authors give for- 
mule for the intensities of the orders 0 and —1. For the intensities of 
the orders +1 and —2, however, they only give expressions in the case of 
exact Bragg incidence. 


The purpose of the present paper is to deduce these results from the 
generalised theory of Raman and Nath? and to give also expressions for the 
intensities of the orders +1 and —2 in the case that the angle of incidence 
is close to the Bragg angle. 


2. THE RESULTS OF AGGARWAL® 


According to the generalised theory of Raman and Nath, the amplitudes 
of the diffracted light waves satisfy the following system of difference— 
differential equations, 


2 EP — ba — bss) = (np — 2an sin P) $n 
(n = — co, ..., +00) (1) 
with the boundary conditions 
$o60)=1, ¢2()=0 (#0) (2) 
where 


__ 2m4Z 
g = 2nd 


§ At present at Rijksuniversiteit te Gent, Laboratorium voor Kristalkunde, Ghent, Belgium . 
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- & 
P= sek 
— A 
a= a 


A = the wave-length of the incident light in vacuum, 

A* = the wave-length of the supersonic waves in the considered medium, 
fy = the refractive index of the undisturbed medium, 

#, = the maximum variation of the refractive index, 

Z =the distance travelled by the light beam through the ultrasonic 


field, 
y = the angle of incidence of the light. 
Supposing $;. = ¢,3; = --- = 0, Aggarwal has shown that if it is assumed, 


that the intensity 7, = 1 and remains so when the lines of order + 1 and — 1 
appear, the expressions for the first order intensities are given by 


me (my sin? Ee (ae + 2 sin )| 
with s., E ( + 2 sin v)| 
2A* \uA* 
L being the width of the sound field.+ 


We notice that the expressions (3) remain finite when the angle of incidence 
is equal to the Bragg angle + for the order — 1, defined by 





(3) 


a 
sin fy = 4 aw? = 


, AL 
Iy=L, amt (Fo) 
(a) 


One gets 


Now the intensity of order + 1 may be neglected if aan ak. &e & & 


* The expression for p in the Raman-Nath theory contains a factor cos?%. Resuming the 
reasoning of those authors to obtain the system of difference—differential equations, but expressing 
the periodicity along the direction of propagation of the sound waves instead of along a direction 
normal to the incident light beam, one obtains the above value of p. 

+ In fact we should have to read L/cos% instead of L in Formula (3). In practice, however, 
the angle » is always very small, so that the distance traversed by the light in the ultrasonic field 
may be approximated by L. 
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evident that while for low frequencies of the ultrasonic waves the orders ++ 1 
and — | appear simultaneously, for higher frequencies only the order — 1 
appears. J, remaining equal to unity and /_, becoming maximum for Bragg 


2 
incidence, i.e., I_, =(%) is not very small, it is clear that when the angle 


of incidence is close to the Bragg angle, the expressions (3) are not more valid. 


3. INTENSITIES OF THE ORDERS 0 AND — | 
From the preceding section it is evident that we should obtain better 
expressions for the intensities of the orders 0 and — 1, in the case of an angle 
of incidence close to the Bragg angle, if we assume that in the system (1) all 
the amplitudes, except ¢, and ¢,, disappear. The system now becomes 


7 dby 


+ ¢,=0 
2 (4) 
2 Oe — by = ding, 
where 
4o = p — 2a sin & (5) 


the boundary conditions still are 


do (0) = 1, $y (0) =0 | (6) 
The equations (4) form an easy system of differential equations of the second 
order, the solution of which, satisfying the boundary conditions (6) is, 


by = [cos Vor 4: E—- Ter sin Vt FH- elect (7) 
Bi sini nid vets x J eice (8) 


The intensities of the orders 0 and — 1 immediately foilow from (7) and (8). 
After effectuating the calculations we get,{ 


ly = deta = ary [o* + $008? VF FE: €) ) 
Ly, = $a = 4 SEE (10) 


t Generally one has for the intensity of the order n, J, = ¢-» ¢-n* in accordance with the 
preliminary theory of Raman and Nath, where /, is found as the square of a Bessel function of 
the order — n. 
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4. 


In order to compute this intensity, we shall follow Aggarwal’s method: 
we assume that the expressions for the intensities of the orders 0 and — 1 
do not change when the line of order — 2 appears. The amplitude of the 
order — 2 can then be determined by the following differential equation of 
the system (1), 


INTENSITY OF THE ORDER — 2 


d ’ : 
2 2 — $, = i (4p — 4a sin #) $4 (11) 
The solution of this differential equation satisfying the boundary conditions is 
eke eke 
rs hee 4 | xk (k,—ks) + (k,—ke) (ks—kg) 


etk.g 
+ ee] ™ 


where 

k, =cCo oa Vo + 4 

kg =o — Vor?+4 

k, = 2(p — asin #) (13) 
The intensity of the order — 2 consequently is given by 


12 = doh2* = [Ke — ks) (ky — Ka) (ky — Ke)" 


a (ke—ke) : sin? (ks—k,) : 
k, _— Ks ” ks tie ky 








sin? (k,—k,) : 
+ pee (14) 


Now if we only consider an exact Bragg incidence, so that 
o=0, ki =}, kz = — 4, kg =p, 
the formula (14) becomes 
! sin*(p+4)5 sin*(p—p§ | 
Aon pal “Ste ee Pee 
It is well known‘ that if only a few order lines appear in the diffraction spec- 


trum, p>>1; so 4may be neglected with respect to p and the formula (15) 
finally becomes 












It 


If 


Tal 
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ine o sin? § (16) 


5. INTENSITY OF THE ORDER + 1 


In the same way as in Section 4, we see that the amplitude of the order 
+ 1 is determined by the differential equation 


22 + dy = ilp + asin Ye, (17) 
Its solution satisfying the boundary conditions is 


Beenie (eh? 4 __ kena 
he 2 L(ki—Kk4) (ka—ky) © (Ka—ky) (Kg—k) 


(18) 





keke 
+ Rk) (kk) 


where k, and k, have the same values as given by (13) and 

k,=4(e + 2a sin #) (19) 
Consequently the intensity of the order + 1 is 

Ty = $4 $4* = 4 [Ki — he) (i — fea) (Ka — I? 


= sin? (kk) § kg (Kat hak.) sin®(ky—ky) § 








k, — ky kK, — ky 





Ak, (ky+a—K,) sin? (ky Ka) | 
* ka — Ke 
(20) 
If again we consider an exact Bragg incidence, so that 


o=0, k, = 3, k, = — }, k,=p 


this rather complicated expression becomes 





Ia = 4° iy {— sin + 72, sint(o— 9 § 


2p gine £ 
+ Py sinte+H5} 20) 
Taking into account that p> 1, this formula reduces itself to 


hh: ss [ — sin® 5+ 2.5 sin (p + ¥ 5 | (22) 


A2 
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or 





l — . 5 pe 
l= 4-5 [- sin? 5 + 4 sin® 5 | (23) | 
where the symbol ¥ has the following meaning, 5 f(a+b) = f(a+b) 
+,- +,- 
+ f(a— 5). 
6. COMPARISON WITH THE RESULTS OF BHATIA AND NOBLE 
There exists a simple relation between the maximum variation of the 
refractive index », and the maximum variation of the scattering index 4 
used by Bhatia and Noble, 
4 = #01, where a = $(j49? + 2) (uo? — 1) = 1 
i 
Using this relation and introducing the notations of Bhatia and Noble, : 
( 
vo A*\? _ ma __ pe? .. | 
8=4(5). B= ye f="y Sin ¢ : 
(it has to be noticed that this has another meaning than the same notation te 
used above) in the expressions (9), (10), (16) and (22) we obtain exactly li 
formule (4.12 a), (4.12 5), (4.12 c) and (4.12 d) of Part II of their paper. th 
It must be noted that formula (4.12 d) of those authors contains two print- ct 
ing errors, which were kindly confirmed to us by Dr. A. B. Bhatia. a 
7. SUMMARY : 
It has been shown that in the case of oblique incidence of the light when 
the angle of incidence is close to the Bragg angle and for high frequencies d 
of the supersonic waves, the intensities of the orders 0 and — 1 and of the | 
orders + 1 and — 2 can easily be obtained from the generalised theory of 
Raman and Nath. The formule are the same as those obtained by a more SI 
involved perturbation calculation by Bhatia and Noble. Sp 
We are very thankful to Dr. R. Mertens for the many fruitful discussions th 
we had with him on this subject. by 
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INTRODUCTION 


1. In this paper I prove theorems relating to the products of matrices 
which define either convergence preserving or limit preserving transforma- 
tions of a series or a sequence into either a series or a sequence. The pro- 
perties of sequence-to-sequence, series-to-sequence and series-to-series trans- 
formation matrices have been studied earlier by several authors and a com- 
prehensive account of them can be found in Dienes*? and Cooke.! The 
sequence-to-series transformation has been defined by Hill.* The products, 
1 taken two by two, of some of the above matrices have been studied by Vermes.*-? 
y In Sections 3 and 4 of this paper, I have considered the remaining cases of 
the products. As a consequence of this, I prove that the matrices defining 
convergence preserving sequence-to-series transformations form a Banach 
algebra under a suitable norm. Also, in Section 5, I prove two simple 
sufficient conditions for a suitably restricted sequence-to-series method to 
include or be equivalent with another of a similar type. 





: 2. Notations and known theorems.—I shall briefly recall the various 
“ definitions and theorems proved on the products of these matrices. 

of We shall denote by (c) the space of convergent sequences, by (cy) the 
e space of null sequences, by (y) the space of convergent series and by (y,) the 


space of convergent series with sum zero. The matrix A= ay, transforming 
the ‘ points’ of the space (X) into those of (Y) shall be, for brevity, denoted 
by A c I'(X, Y). 


The standard method of summing a series uy + uy, + Ug +.......... 
with partial sums {s,}, sp = Up + uy+....-+-Uz is the transformation of the 
sequence {s;,} by the matrix A= ay, into a convergent sequence on= J ankSx. 

k 


The matrix A is said to be convergence preserving or simply a K-matrix if 
the convergence of {s;,} implies that of {c,}. If in addition the limit is also 
preserved then the transformation is said to be regular and the matrix a 
T-matrix. The matrix is said to be “ multiplicative p”’ or simply c I (c, pc) 
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if lim s; = / implies lim on, = p./. The necessary and sufficient conditions 
that the matrix A is a K-matrix are (Dienes*, p. 388) 


(2.1.1) 3 | anx | <M, a quantity independent of n, 
k 


(2.1.2) limn ank = 5, exists for each k, 
and (2.1.3) yn= JZ ank—y as n — 00. 
k 


ee 


The necessary and sufficient conditions that A is “ multiplicative p” are 
that 5; =0 and y =p; consequently A is a 7-matrix if and only if, 5, =0, 
y= 1. Conditions (2.1.1) and (2.1.2) above are necessary and sufficient 
for A to be ¢ I'(éo, c) while (2.1.1) and (2.1.2) with 5, = 0 are necessary 
and sufficient for A to be C I'(co, co). The matrix A satisfying condition 
(2.1.1) shall be called a K,-matrix. 


Another method is the transformation of the series Su, into a convergent 
k 


sequence tn = J fnkuk by a matrix F=fpx. A convergence preserving 
k 


series-to-sequence transformation matrix is called a B-matrix while the limit 
preserving one is called a y-matrix. The necessary and sufficient conditions 
that F may be a f-matrix are: 


(2.2.1) 3 \fnk —Sn, kia | <M, @ quantity independent of n. 
k 


and (2.2.2) limnfnk = Bx exists for each k. 


The matrix F will be a y-matrix if, and only if, in addition £;=1 and F 
will be c r(y, Co) if B, = 0. 


A third method is the transformation of the series 5 uj, by the matrix 
k 
H= hy into a convergence series S vpn, so that v;2= J hyyu,. A convergence 
k 


preserving transformation of this type is said to be defined by a 8-matrix 
and a limit preserving one by an a-matrix. The necessary and sufficient 
condition that the matrix H is a 5 or an a-matrix is that the matrix F defined 
by 


(2.3) fnk = how + Puk +----+ Ank OF hank = Snk — fra, k 
is a B-matrix or a y-matrix, respectively (Vermes® ”). 


The sequence-to-series transformation, defined by Hill,? consists in 
transforming a sequence {s,} into a convergent series SY vp, by a matrix P= pax 
n 


so that u%, = Y PnesS~. We define a matrix yielding a convergence preserving 
; k 


- 

















it 
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transformation of this type to be a A-matrix while the limit preserving one 
is called a p-matrix. The necessary and sufficient condition for the matrix P 
to be a A-matrix or a p-matrix is that the matrix A defined by 


(2.4) Gnk = Pok + Pik +----+ Pnk OF Pnk = 4nk — In, k 
is a K-matrix or a 7-matrix, respectively (Hill*). 


From the definitions and properties mentioned above, we have also 
THEOREM 2.5: Every 5-matrix is also a B-matrix and indeed < I'(y, co). 
Also, in a A-matrix P, we have, if A is the corresponding K-matrix then 


- | Pnk | =2 | Qnk—4n-1, k | <2 | nk |+ 2 | ana, i | <2M, 


limp Pnk = 0, since limy an exists. 
Also, 
=z Pnk = 2 (Ank — 4-1, k) = Yn — Yn-a 
k k 


and therefore 
limy, ~ Puk = 9. 


Thus we have 

THEOREM 2.6: Every A-matrix is a multiplicative zero matrix; i.e., 
eI'le; ea). 

Several theorems on the products of the matrices have been proved earlier 
and we shall briefly state them here for the sake of completeness. 

THEOREM 2.7: The product C = AB of two K-matrices is a K-matrix. 


THEOREM 2.8: The product C = AF of a B-matrix F with a matrix A 
exists and is a B-matrix if, and only if, A is a K-matrix. 

THEOREM 2.9: The product C = AF of a y-matrix F with a matrix A 
exists and is a y-matrix if, and only if, A is a T-matrix. 


THEOREM 2.10: The product C = FA of a y-matrix F and a T-matrix A 
need not exist. 


THEOREM 2.11: The product C = FH of any B-matrix F with a matrix 
H exists and is a B-matrix if, and only if, H is a 5-matrix. 


THEOREM 2.12. The product C = FH of any y-matrix F with a matrix 
H exists and is a y-matrix if, and only if, H is an a-matrix. 


THEOREM 2.13: The product C = GF of two B-matrices G and F need 
not exist. 
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THEOREM 2.14: The product C = HF of an a-matrix H with a y-matrix 
F need not exist. 


THEOREM 2.15: The product of two 5-matrices is a 5-matrix. 
THEOREM 2.16: The product of two a-matrices is an a-matrix. 
THEOREM 2.7 is well known and its proof can be found in Cooke. 
THEOREMS 2.8-2.16 are due to Vermes*’. 


3. In this section I shall prove some preliminary theorems on the pro. 
ducts of summability matrices. 


We have in virtue of the various necessary and sufficient conditions 


stated earlier, the following relationships between the various summability 
matrices. 


BCACT(e,¢)) CK CB andaci cI(y,¢9) cB. 


Also, from the way the A and 5-matrices are formed out of their corresponding 


K and f-matrices we have, from the relation K c BthatA cd. Buty ca 
since T ¢ y. 


THEOREM 3.1: Let A c I'(c, pc) and B c I'(c, pc). Then the product 
C = ABis c I'(c, pp’c). 


Proof.—It follows from the hypothesis that, since both A and B are 
K,-matrices, so is their product C. 


Now, ¢nk = z Anibix, and therefore for each fixed k, {cy} can be regard- 


ed as the A-transform of the sequence {b;,}. 
But, since {b;;,} is a null sequence, by hypothesis, and also since 
A cI (c, pc) we have that limncn, = 0, for each fixed k. 
Further, 5 cn, = 5 DS anibix- 
k k i 
Since A and B are K,-matrices, we have that Y ¥ | ani|-| dix | < 00. 
k 4 
Therefore we have, on the interchange of the order of summation, 
=X Cnk = FF ani dix = Z Ani - biz = 2 Aniyi, Say. 
k 4 ek i i 
But since B c I'(c, p’c), yi’ —>p’ as i—> 00. 
Zcnk can now be regarded as the A-transform of the sequence (y;’) and conse- 
k 


quently we have that limn Y cnk = p-p’- 
k 


This completes the proof of the theorem. 




























THEOREM 3.2: 
H is a B-matrix. 





so that 


need not exist. 


previous theorem. 
THEOREM 3.4: 





The product C = AH of a K-matrix A and a 8-matrix 





fo Seer 
iE eer 
1 ae 
2 8 
Sai 
AH = 





The product C = AH of the matrix A < I'(c, pc) and 
d-matrix H is < I'(y, Co). 


and H = 








coe eeee® 


which is neither a K- nor a 5-matrix. 
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The conclusion of the theorem is seen to be true immediately from 
Theorem 2.6, since every 5-matrix is also a f-matrix, by theorem 2.5. 
Now, we shall prove that the conclusion is the “‘ best possible”, in the sense 
that it need be only a f and not necessarily a K- or a 5-matrix, since it is likely 


that a B-matrix may be a K- or a 5-matrix. This is seen by taking 


eee eee eee eeee 


Note.—For the product C to be a f-matrix for every K-matrix A, 
it is not necessary that H should be a 5-matrix, since by Theorem 2.8 it is 
enough if it is a B-matrix. 


THEOREM 3.3: The product C = HA of a 5-matrix H with a K-matrix A 


This can be directly verified by taking the matrices A and H as in the 


Proof.—Since any multiplicative matrix is also a K-matrix it follows 
from Theorem 3.2 that the product matrix C is at least a f-matrix. 
it is enough if we prove that limncn, = 0. This is evident since, as before, 
we can regard the sequence {Cn} for each fixed k, as the A-transform of the 
sequence {hj,}, regarded as a sequence in /, which in the present case happens 
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Now, 





176 M. S. RAMANUJAN 


to be a null sequence on account of the convergence of {nk = Aok+ Ay, + 
.... + hae as nm —>co and also since the matrix A is multiplicative. 


COROLLARY: The product C = AH of a T-matrix A and an a-matrix 
H is cI (y, Co). 


That the product need not be a 5-matrix can be seen by taking for A 
the (C, 1) matrix and H, the Identity Matrix, so that the product is still the 


(C, 1) matrix, which is not a 5-matrix. Thus the conclusion is the “ best 
possible”. 


4. A and p-matrices.—In this section we shall consider the products 
of the various matrices with the A or w-matrices, either to the left or to the 


right. 


THEOREM 4.1: The product C = PA of a dA-matrix P with a matrix A 
exists and is a A-matrix if, and only if, A is a K-matrix. 


Proof.—Let B be the K-matrix corresponding to the A-matrix P, 
so that bnk = Pok + Pik +----+ Pnk OT Pak = bnk — Dna, k- 
Now Cnk = 2 Pni Uk = 2 (oni — Ona, i) dik = [BA]nk — [BA}n-1, k: 
But since B and A are K-matrices, so is their product BA, and therefore it 


follows that C is a A-matrix. Thus the condition is sufficient. That the 
condition is necessary is seen by taking 








: 8° Bios 
—2 2 © ‘Bhismssus 
O82 & Girciacds 
P= and any matrix A. 
S BB Boivcscns 


Then Cnk = ank — Aan-1, k- 
Now, if C is to be a A-matrix, then A must be a K-matrix. 


The above proof also gives rise to 


THEOREM 4.2: The product C = PA of apu-matrix P with any matrix A 
exists and is a y-matrix if, and only if, A is a T-matrix. ' 


THEOREM 4.3: 
P is a K-matrix. 


The product C = AP of any K-matrix A with a d-matrix 











e 
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The result follows from the fact that every A-matrix is also a K-matrix and 
Theorem 2.7. 


We shall also prove that the above conclusion is the ‘ best possible’ 
by proving that the product need not be a A-matrix. This will follow by 
taking for A the matrix A of Theorem 3.2 and P any arbitrary A-matrix. 
Then, 

Poo Poi: -Pok------ 


Poo Po: -Pok------ 


eoeoreerereeereeeereeese 








and this cannot be a A-matrix, unless {py;} are all zero. 

Remark 1.—The product C = AP of the matrix A c ['(c, pc) anda 
\-matrix P is c I'(c, ¢9). 

This is a consequence of the fact that every A-matrix is I"(c, co) and 
Theorem 3.1. 

From Remark 1, it follows that 

Remark 2.—The product C = AP of a T-matrix A and a w-matrix P 
is ¢ I'(c, co) and therefore cannot be a T-matrix. 


Remark 3.—In order that the product C = AP of any K-matrix A with 
a matrix P may be a K-matrix it is not necessary that P should be a A-matrix, 
since by Theorem 2.7 it is enough if it is a K-matrix. 


Remark 4.—In order that the product C = AP may be a K-matrix for 
any arbitrary A-matrix P and a matrix A it is not necessary that A should 
be a K-matrix as may be seen by taking for A the matrix 








eoeeee eee eee ees 


Then Cnk = J 4niPik = Zz Pik=5x Say, where 5; is the quantity lim; bj, when 
> 


B is the K-matrix corresponding to the A-matrix P. Now the result follows 
from the fact that for any K-matrix Y | 54 | < oo. 
k 
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THEOREM 4.4: The product of two d-matrices P and Q is also a A-matrix. 


Proof.—Let A be the K-matrix corresponding to the A-matrix P defined 
by the equation (2.4). Then, we have, if C = PQ 


Cnk = 2 Pnifik = F (Ani — An, i) Fik = [AQInk — [AQIna, k 


and AQ is by Theorem 4.3 a K-matrix. Therefore the result follows. 
From Remark 2 to Theorem 4.3 we have, 


THEOREM 4.5: The product of two pu-matrices P and Q is c I'(c, y) 
and cannot be a w-matrix. 


THEOREM 4.6: The product C = PF of a d-matrix P with amatrix F 
is a 8-matrix if, and only if, F is a B-matrix. 

The result is a consequence of the fact that if A is the K-matrix corres- 
ponding to the A-matrix P and D = AF then dyp=—cont+ Cixet+.----+ Cnk 
Now C is a 5-matrix if, and only if, D is a §-matrix, which happens by 
Theorem 2.8 if and only if F is a B-matrix, since by hypothesis A is a K-matrix. 

A similar argument shows that 

THEOREM 4.7: The product C = PF of a p-matrix P with a matrix F 
is an a-matrix if, and only if, F is a y-matrix. 

THEOREM 4.8: The product C = PH of a d-matrix P with a 8-matrix 
H is a 8-matrix. 

Proof.—Let A be the K-matrix corresponding to the A-matrix P. Then 
as before, Cn~ = [AH]nk — [AH]n-1, k- 


Now, AH is by Theorem 3.2, a f-matrix and therefore we have that C is a 
5-matrix. 


It may be remarked here also that in order that the product C = PH 
of a matrix P with a 6-matrix H may be a 5-matrix it is not necessary that P 
should be a A-matrix, since we have by Theorem 2.15 that it is enough if 
it is a 5-matrix. 


After the Corollary to Theorem 3.4 we have also the following: 


THEOREM 4.9: The product C = PH of any p-matrix P and an a-matrix 
H is Cc r(y, Yo)- 


Next we prove the following theorem. 


THEOREM 4.10: The product C = FP of a B-matrix F with a d-matrix 
P is a matrix c I'(¢9, c). 
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Proof.—Since every A-matrix is also a K-matrix, we have, after Theorem 
2.11 that C is a B-matrix; i.e., that each column of the matrix C is con- 
vergent. 


Now if A is the K-matrix corresponding to the A-matrix P, we have 

Cnk = 2 SniPix 

= LSni Gik — Gia,n) = & Sri — fn, ivr) Gi + lim; (Snidix)- 

Therefore, we have 

- | Cnk |< = | 2 (Sni — fn, isd dip | + x | lim; ( fné ix) | 
But since F is a B-matrix, we have that {/,;} are bounded, regarded as 
sequences in n or i (Vermes’). 
Therefore, 

= | enk | <22 \fni —Snjinr | lai | + =M | 8x | 

k i 
where M is the upper bound of the sequence {fi} regarded as a sequence 
in i and lim; diz = 5x. 
Since now in a K-matrix ¥ | 5, |< oo we have that J | cn_, |<oo, for each 

k 


n, since F is a B-matrix and A is a K-matrix. 


Thus the product C is also a K,-matrix. Now the conclusion in the theorem 
follows from C being simultaneously 6 and K;. 

THEOREM 4.11: The product C = FP of a y-matrix F and a w-matrix P 
is a matrix C I’ (Cp, Co). 


Proof.—In the light of the proof of Theorem 4.10, we have only to show 
that limncn,= 90. This is obvious since {cp,} can be regarded as the F-trans- 
form of the series 2 Pik which by the property of the u-matrix is a series with 


sum zero. 
Theorems 4.10 and 4.11 give rise to the following two theorems: 


THEOREM 4.12: The product C = HP of a 8-matrix H with a d-matrix 
P is a matrix © I'(€o, y). 


THEOREM 4.13: The product C = HP of an a-matrix H with a w-matrix 
P is a matrix <I (Co, yo). 


THEOREM 4.14: The set of A-matrices P form a Banach Algebra under 
the norm || P || = 2.1.u.b. 5 | S pix | 


k=0 i=0 








180 M. S. RAMANUJAN 


Proof.—Let P” = [pix] be a sequence of A-matrices. 
Then P”.P* and P” + P* are A-matrices is evident. 
Therefore we have only to prove that the space of the above matrices is 
complete under the norm specified and that || P™. PS || < ||P” ||.|| PS\|. First 
we prove the completeness of the space, i.e., that each Cauchy sequence in 


the space implies the existence of a member of the space to which that sequence 
will converge. 


Now, let || P™ — P*|| < «, r, s> np. 


That the above condition is necessary for the convergence of {P"} is clear. To 


prove that it is also sufficient, let A” denote the K-matrix corresponding to 
the A-matrix P’. 


Now, given a matrix D = [dix], let || D ||, denote 2./.u.b 5 | 3 dix | 
n k 


t=0 


and let || D ||, denote /.u.b X | dnx |. 


Then we have that || P* ||, = 2 || AT |l, 
Now, || P™’ — P§ ||, = 2 || A” — A® lk <¢, 7, 52> Mg. 


Since || A7 — A® |i, </2, r,s=>m, and the K-matrices form a Banach 
Algebra under the norm |j A ||,* we have that there exists a K-matrix A to 
which the sequence of matrices A” converge. Corresponding to this K- 
matrix A there exists a A-matrix P which is the limit of the sequence P’. Thus 
the space is complete under the norm chosen and indeed complete under 


the norm || P || = /.u.b 2 | 2 Pit | 

To prove that || P”.P® ||, <|| P” |, . || PS lh, 

lt C=P".P*, i.c., ca, = Z P' ni - Pri 

Let B= (bnx) where bn, = Con + Cry + .--- + Cnk 
ie. B= A’ .P*. 


Now, since A” is a K-matrix and the matrix P* is a A-matrix and therefore 
necessarily a K-matrix we have that B is a K-matrix and therefore 

l| B ilk <|| AT Ile - I] PS ilk 
and consequently || C ||, = 2 || B |lk <2 || A” lx - || P® ll = |] P* Il, - | PS Ile 


* This result is due to Parameswaran‘ (Theorem J). 
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Now, || P® lle = /.u.b 2 | Penk | 
= l.u.b z | ante — na, k | 
Ss l.u.b = | ane | + tub ~ | na, | = 2 || AP lle = |] PS I. 


Therefore ||C ||, <|| P” ll, . || P* ll 
and this completes the proof of the theorem. 
Remark.—On account of the 1-1 correspondence between the A-matrices 
P’ and the K-matrices A” defined by the relation 
nk = Poh + Puke + ---- + P'nk 
and the fact that the K-matrices form a Banach Algebra under the norm || A ||, 
it is natural to expect that the A-matrices will form a Banach Algebra under 


the corresponding K-norm. But that they do not do so can be seen from 
the example below. Suppose we have || P ||, = || A Il, 








Let iO Ow... 
ae ft Pes 
pr = Ps — 0-1 1 O....} so that AT = A’ =I, identity 
matrix 
G: @—}i boc 


In the notation of the proof of the theorem we have, D = P*. 
Therefore || D ||; = || P* ||, = 2. 


But || A” ||; = 1 and therefore || D ||, > || A” |lx . || A® Ilk, i-e., 
WC n> U1 PP ll, - PF th. 


5. Inclusion Theorems.—In this section we obtain a necessary and suffi- 
cient condition for a sequence-to-series transformation method defined by a 
suitably restricted »-matrix P to include another similar method, given by 
a matrix Q. Also a sufficient condition for the inclusion of one method in 
the other is proved. 


Definitions.—Let {s,} be any sequence, convergent or not, and let P and 
Q be two u-matrices. If the existence of the P-limit of the sequence {s,} 
implies the existence and its equality with the Q-limit of the sequence then 
the method Q is said to include the method P and this is denoted by Q > P. 
If in addition P > Q, then the methods are said to be equivalent and this is 
denoted by P~ Q. 
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THEOREM 5.1: Let P and Q be two row-finite u-matrices and let P have 
a left-hand row-finite reciprocal “P. Then Q> P if, and only if, “P isa 
y-matrix. 


Proof.—Let {s;,} be any sequence. 
Let vz = 2 Pnk Sk and XS’ vn = I, say. 
Now, since P has a left-hand row-finite reciprocal -'!P, we have 
Sm = 2 “lpn Um: 
If now Y yu is the Q-transform of the sequence {s;,) we have, 
ui = 2 dim $m = 2 dim Z —Pmn Yn: 
Since 2 vy, converges and Q and —'P are row-finite, we have on the interchange 


of the order of summation of the above expression that the series Yu; is the 
i 


(Q.-'P) transform of the series X' vy. Therefore, if now, we must have that 
Q oP, ie., Xu; = YF vy = / then the matrix (Q.-'P) must be an a-matrix. 
But since Q is a w-matrix, it follows from Theorem 4.7 that Q.-!P will be 
an a-matrix if, and only if, —-'P is a y-matrix and this completes the proof 
of the theorem. 


COROLLARY: Let the u-matrices P and Q with their left-hand reciprocals 
-1P and Q (supposed to be existent) be row-finite. Then P ~ Q if, and only 
if both-*P and~Q are y-matrices. 


Even if we do not assume the existence of the reciprocals, we can still 
proceed as follows to find a sufficient condition for Q> P. Let {s,} be any 
sequence and v, = Y' PnkS~. Let Y vn = I, say. 

: 


Now take the transformation ¥ um of the series Xv, by an a-matrix, H, say. 
Then Um = J hmn Vn = J Amn J Pnksk- 

n n k 
If both H and P are row-finite then the interchange of the order of the above 


summation will be valid and we will have that Xu, is the transform of the 
sequence {s;,} by the matrix HP. 


Therefore, if we have that Q = HP then we have that Q > P. 
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Thus we have the 


THEOREM 5.2: Jf P and Q are two row-finite y-matrices and there exists 
a row-finite a-matrix H such that Q = HP then Q 3 P. 


COROLLARY: Jf P and Q are two row-finite -matrices and there exist 
two row-finite a-matrices H and H' such that Q = HP and P = H'Q then 
P~@Q. 


I would like to draw the attention of the reader to the table I have given 
in the Appendix, which gives a complete picture of the existence and the 
classification of the various products considered in this paper and by the 
previous authors as well. 


In conclusion I have very great pleasure in recording my indebtedness 
to Professor C. T. Rajagopal and Professor V. Ganapathy Iyer for their help 
and guidance in the preparation of this paper. 


Note.—I learn from Dr. P. Vermes of the University of London that 
Mr. F. Gerrish of that University has obtained independently results similar 


















































to those in this paper. I am grateful to Dr. Vermes for drawing my 
attention to the work of Mr. Gerrish. 
APPENDIX 
‘oes | S| Le | 
| Fal B § A | | y ee OT, 
| ele ts K | | tT | » (PO, CoH, Co) 
K | 
01) [5]| (3-2) (43) [1] [5]} (3+4)| (48) 
need not|need not| B P'(Co, C)| need notineed not} y (Co, Co) 
We. exist exist | | exist exist | 
[5] [6]| [7]| (4+10)| | [5] [6] [6]) (4-11) 
| need not|need not| 8 I'(Co, 7) | Ineed notijneed not| a I'(Co, Yo) 
FY exist exist | | | exist exist 
(3-3) (6)| [7] (4-12)| | (3°3) (6) [6}]| (4-13) 
| d &§ | 8 A | Ce a Cy, vo)\I"(C, vo) 
r } | | 
(4:1)! (4+6)| (4-8) (4+4)| : | (4-2)1 (4-7)1 (49) (4-5) 























The above table gives the complete picture of the existence and the classi- 
fication of the various products considered in this paper and also by the 
previous authors. The numbers in the brackets ( ) refer to the theorem 
in this paper where the result is proved, while the numbers in the brackets 
[ ] refer to the literature cited at the end of the paper. The product of any 
two matrices can be found out immediately from the table above. 





184 


1. Cooke, R. G. 
Dienes, P. 
Hill, J. D. 


4. Parameswaran, M. R. 


5. Vermes, P. 


M. S. RAMANUJAN 


REFERENCES 


Infinite Matrices and Sequence Spaces, Macmillan, 1950. 

Taylor Series, Oxford, 1931. 

**On the Space (y) of Convergent Series,” Tohoku Math. Jour., 
1939, 45, 332—37. 

**Some Converse Theorems on Summability,” Proc. Ind. Acad. 
Sci., 1952, 36, 363-69. 

“Product of a T-matrix and a y-matrix,” Jour. Lond. Math. 
Soc., 1946, 21, 129-34. 


“* Series-to-Series Transformation and Analytic Continuation by 
Matrix Methods,” Amer. Jour. Math., 1949, 71, 541-62. 

** Conservative Series-to-Series Transformation Matrices,’ Acta 
Scientarum Mathematicarum (Szegad), 1951, 14, 23-38. 








—_ > 2 soe 














HYPO-ELASTIC PURE FLEXURE 


By P. D. S. VERMA 


(Department of Mathematics, Indian Institute of Technology, Khargpur) 
Received May 31, 1956 


(Communicated by Dr. B. R. Seth, F.A.sc.) 
1. INTRODUCTION 


IN hypo-elasticity, a subject of very recent origin, Truesdell’ has generalised 
the classical linear theory of elasticity by taking the stress increment as a 
function of small strain from the immediately preceding state. Mathema- 
tically speaking his constitutive equation is 


Rate of stress = F (Rate of deformation). (1.1) 


For the rate of stress he takes, to be dynamically consistent, the tensor 
~ > ? : , : 3 
Sj} = 5, Sib + Si'o®, ie — S05, 1 — Si*0*, ie + Syv® (1-2) 


where S;* is the spatial stress tensor referred for convenience to twice the 
shear modulus as the unit of stress, v', the velocity, comma, the covariant 
derivative; and for the rate of deformation 


dij = 4 (ri, 5 + %, i). 
In a subsequent communication, (Truesdell’), he takes a still more general 
concept of elasticity by writing the constitutive equation of the ideal material 


as S = F(d, S), where response to deformation d is moderated by_ reigning 
stress field S. Assuming that no modulus of the ideal material shall carry 
a dimension independent of the dimension of stress and that F is an ana- 
lytic function of the components d;‘ and S;', he establishes the following 
differential equations for an hypo-elastic isotropic body: 


S = 8gol + gd + 59,8 + Mggl + 4g, (dS + Sd) 
+ 8g,S* + Mg.S + Ng7l + 4g, (dS? + Sd) 
+ Mg,S* + NgioS + NgyS?’, 


where g’s being dimensionless analytic functions of the principal invariants 
of S, and for different grades, g’s being defined in a particular problem. 
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2. Here three equations of motion and one equation of mass have 
to be taken as in the classical theory of elasticity. Boundary conditions are 
also the same. The most important and interesting feature of the above 
non-linear theory is that the stress-strain relations are not assumed. 
Truesdell? has already discussed problems on simple shear and torsion. 


3. In this paper, we propose to discuss the steady bending of a narrow 
but thick rectangular plate of constant density into a region of circular form 
of inner radius ‘a’ and outer radius ‘ >” assuming the material to be hypo- 
elastic and isotropic of grade zero. When the compressible material reduces 
to an incompressible one, it is interesting to find that the solution is inde- 
pendent of the assumption of the stress-strain relations. The results obtained 
here are found to be identical with those obtained by Seth® in which the 
assumption of stress-strain relations is involved. 


4. For an isotropic hypo-elastic body of grade zero, we have: 
A 
80 = Fy 8, =1, 28 =8:=8.=....-.. = 21, = 0, for all S. 


In cylindrical system of co-ordinates, the six constitutive differential equa- 


tions of hypo-elasticity for a body of grade zero and of constant density p 
reduce to 


Py) ry d 0 du 
[st ap +? 50 + ™ 55 | Se” — 28" 5 


or 
-asyMosy(1+Ay= Bae 


=4[QteSh @2) 


Q Q Oo M1or_g9M@_o2 MU 
[setae t? gat ” 55] 8 Sz 00 Sz dz 


ow ou 1 dw ow 
= T it = _— a Se _ ieeseiend 
Sz az T >) So r2 00 Sr ar 


=1(S+ 7} 49 











La + "os 


Ls + 





” 
2u [pS + 5 30 Se" + 


So + 5S" + - 


2h E 
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om ytws 7 | Se — So Sr? +1) 
99 (2+ 3)-searen a(t 8) 0 
ee ae | Se? — S2(; dz sta 


1 dw 
ae : ls = “VWaw-* >; 


-4[Z+ a} «9 


3 ry 3d Lh zw r OW 
us, + 39 +5. |S? 28? sz — 282" 5, 


— 5F o ( ? n) =; a itt 


The three equations of motion in the absence of external forces reduce to 





S,” ea S,° ri) 
ar CS S." | 


=p|u™ -+ ¥ g-) two sh (4.7) 


=p[us +e Se a ic (4.8) 
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In pure flexure, the shears S,” = S,’ = S,? =0, and we tentatively 
take S,”, S,°, S as functions of r only. Considering (4.1), we get: 


d d 
u 5, Se” = (1 + 28,7) 
which gives: 
u= x’ (0) /1 + 2S,", (4.11) 
where x’ (9) is any arbitrary function of 6. Considering (4.4), we have: 
d .9 9. fu , w 
u5Se=(1 + 28 6) (= + 4 
or using (4.11), 
dv l da, 1 
39 = X' (9) VI + 28,7 les ar Se’ — 7} 
or v = x (9) A(r) + F (7), where F (r) is an arbitrary function of r only, and 
—e 1 d 1 
A) = VIF 28" [15559 550? — 3] (4.12) 
From (4.8), we get: 


dv 2uv 
tis cata ‘ats 


Putting for u and v we get: 
7 TT 4e7 » 24 2A(") 7 95e7 
0=x(O[ VIF Isr xX H+AMP+ BO vi +285" | 


+[F@ VI FS" +FOAD +2 viF25,"]. 
For this to hold, we must have: 


2A al 


VI + 28,7 ¥ (rn) + ADP + V1 + 25," =0, 


and 


VIF" F (+ FOAM + =O viFasr =o. (4.13) 


Thus we conclude that F (r) = AA (vr), for all r. 








H 


Sr 
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Hence 
v = [X (6) + AJ A(v). (4.14) 


Considering (4.10), we have: 


' d l J/1 +28," | A(r) 
x’ (#) - | — S,? - ———, 4+ oY “I G4 ee = 
@ La S V1 + 2S," r r | — 


Hence a must be zero. Then either X’ (#) = 0, which gives u = 0; and 
S,"” becomes indeterminate. Therefore we do not take 


x’ (8) = 0. 
Or 
d l /1+ 28,7  A(r) 
So eee Se Beh i as = . 
a” Sita T r + r sa 33) 
Considering (4.2), we get: 
0 Pi) 
- 9 (s,° + +) =r x (4 + S,”), 
or 
— x" (9) {2S.°+ 1} _ _ i), 
%X (8) + A} V1 + 28,0 ~ r2X’ (r), (4.16) 
For this to hold we should have: 
x” (@) = 0, X’ (@) = constant, (4.17) 
r (r) =0. (4.17) 
Using (4.13), we find: 
AM [A® +2 vI+25"] =o. 
Either, 
A=0, 
or 
--2 ym 
A= —< v1 +25; (4.18) 
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If A = 0, from (4.12), we get: 
a. Sep 3 

1+ 28,9 dr a; =O, 

or 
2 

$9 =5 vig ae @, (4.19) 
where D is an arbitrary constant. 
Also from (4.15), we get: 

1 doe, V1+2," 
aces tS 
rV/1+ 2S," =D’, 


= 0, 


where D’ is an arbitrary constant 
or 
D’”? 

Ss," = as 3. (4.20) 

In this case, : 
_ os yee ys 

u= V1 +28," x’ (8) = x’ (0) —, 

where x’ (8) is a constant. 


v =0, and from (4.5), w is also zero. All other equations are satisfied 
except the equation (4.7), which gives 


2p [-?; + a | =>, [ & (0)}? - 4 (- ar) 


this is possible only when D = 0 and then ,’ (@) = a/ 2 (4.21) 
pP 


~ a : VI + 2S,", refer to (4.18), 


then (4-12) gives: 


eS oe Pe 
— 7 Tessar’ — > 
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or 


j,i 
~~ T+ 2860 dr 50% 


or — 2 log r = log (1 + 2S,%) — log A’, where A’ is an arbitrary constant. 


wer 

Se? = 5 — 3 (4.22) 
From (4.15), 

d 1 Vita 2. jee 

le VI+ ost tp OU VI + 2S? = 0, 
or 

d 1 = B”2 

> ll Vit 28 7 V1 + 2S,", or 8,7 = > 7-4, (4-23) 


where B’ is an arbitrary constant. 


From (4.7), we get: 
[Ber fame — SI 
= p [{x' (8)}? rB’? — 2B’ {x (8) + A} x” (6) rB’ 
— r{x (9) + A}? 4 BY], 
which is not satisfied. 
Thus we have the system of stresses and displacements: 


s 7 — D® 1 r) 2 i 
r= a_i So? = — 3,87 =0, Ss = So = S* = 0; 


D’ L 

u= 7 [#0 =0, w=0. 
r p 

Taking the contravariant stress components, we find: 


grr _ DE 4, $00 — 


~ Dr? 


1 
— 77» S* = 0, S76 = S6z — Sz? — (C9. 
Let us apply the boundary condition 
S’7 = 0 at r =a, which gives D’ = a, 
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Thus 


1 


2 
grr — na — 4, $0 = — =e Szz = 0, 


$92 = §'2 = S76 = 0. 


This shows that a traction must be applied to the straight edges and to 
one of the curved surfaces of the cylinder so as to keep the plate in circular 
form. The force which must be applied to the curved surface is 


a 
(S")-» = 49 — 3 


If M, denotes the moment of the couple about an axis perpendicular to the 
plane of the tube and passing through its centre, 


b b 1 b 
M, = f rs?dr= — Jr x adr = — } log ~. 


If — Mg, denotes the moment of the elementary couple in the axial plane 
applied to the plane ends of the cylinder between @ and 6 + d6, we must 
have: 


f rSzdr=0, and — M, =O, 


which is true; 


the traction across the plane 6 = constant is 


1 
0 ax — »- 
S 5,2" 


5. Seth*® has discussed the pure flexure problem in the finite strain 
theory. By taking large strains and assuming the stress-strain relations he 
gets the same stresses as follows: 


2 m 
greater eat Gt! 


9 witp~ O38 


r2 r° 


sre = S’2 = S62 — S2z7 — 0, 


where 





pat = Ot OY (5.1) 
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and 





Gr+y1=C+oe (5.2) 


In the incompressible case, A—> co, c—>0, and using (5.2), 


(4-A+tploA+4yp, ie, 12, 


A? A? 
ST=—ptp og, SV=p—U-g, 


S72 = S79 = Sz = S77 = 0. 


Thus we find that our S’” and his S"” are the same, keeping in view the fact 
that we are taking twice the shear modulus as our unit. S?? is zero in both 
cases. Our S9 and his S® differ only by a constant. We find that the 
traction on one of the curved surfaces and on straight edges must also be 
applied in addition to the couple on the straight edges of the plate whose 
moment is — 4 log (b/a). 


6. By using the strain energy function, Ericksen* has deduced the 
results similar to those given in (4.25). In hypo-elasticity, where we take 
rate of stress as a function or rate of deformation, the stresses are obtained 
without the introduction of a strain-energy function. 


In the end, I must thank Prof. B. R. Seth for his helpful guidance in 
preparation of this work. 


SUMMARY 


In hypo-elasticity, where we take rate of stress as a function of rate of 
deformation, i.e., stress increment as a function of small strain from the 
immediately preceding state, the problems are solved without assuming 
the stress-strain relations. Here we have shown that in the case of bending 
a rectangular plate into circular form, the stresses obtained without assum- 
ing stress-strain relations are the same as those obtained by Seth® after 
assuming the stress-strain relations. 
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ABSTRACT 


The set-up of ‘ Slow-Fast’ conicidence scintillation spectrometer is 
described. Gamma-gamma coincidences in Cs!™ have been studied. The 
604-kev. gamma transition is found to bein cascade with 460 + 20. 
555 + 15, 794 + 15 and 1349 + 30 kev. gamma transitions. The 794-kev. 
gamma transition is found to be in cascade with 604 +. 10 and 555 + 15 
kev. gamma transitions. The results are consistent with the decay scheme 
of Cs!* proposed by Forster et al." 


INTRODUCTION 


COINCIDENCE studies with scintillation counters have been found to be a 
very useful method of determining the sequence of energy levels in a nucleus 
formed in a complex beta-decay of a radio-active nucleus. For this purpose 
a ‘ slow-fast’ coincidence spectrometer as described below has been set up. 
This instrument was used to investigate the gamma transitions in Ba!*4, 
resulting from the complex beta-decay of Cs’. This isotope has been 
studied by several workers’ and in particular, the gamma-ray energies 
were determined from the photo-electron spectrum taken in the Siegbahn- 
Slatis beta-ray spectrometer’ at this Institute. In the decay of Cs!, eleven 
gamma transitions of the following energies were reported—467, 553, 571, 
607, 794, 1027, 1164, 1368 and 1401 kev. Of these the two gamma-rays 
at 604 kev. and 794 kev. are the most intense ones. The aim of this work 
was to determine which of the other gamma transitions are in cascade with 
these two well-known gamma-rays in Ba**, 


In a previous report on beta-gamma and gamma-gamma coincidence 
study of this isotope, G. Bertollini et al/.°, employed discriminators, one 
in each channel, for energy discrimination, the resolving time of the coinci- 
dence unit being 0-5msec. In this work, gamma-gamma coincidences 
have been studied, employing a shorter resolving time, 1-5x10-* sec. of 
the coincidence unit and single channel pulse analysers, one in each channel, 
for energy selection. This helps in deciding more definitely the energies 
of gamma-rays that are in cascade and throws some light on the question 
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of the sequence of levels in Ba?™. In particular, as a result of this investi- 
gation, it can be argued that the first excited state in Ba!™ is 604 kev. rather 
than 794 kev. 

EXPERIMENTAL ARRANGEMENT 


The block diagram of the conventional ‘slow-fast’ coincidence 
scintillation spectrometer is shown in Fig. 1. For gamma-gamma coinci- 
dence measurements two Nal (TI) crystals each 1” in diameter and 1” thick, 
viewed by two EMI6260 photomultipliers have been used. Both the 
crystals were shielded with 4” thick lead collimators. They were kept 120° 
apart so that they do not see each other thus minimizing the false coinci- 
dences arising from scattering of gamma rays from one crystal to the other. 
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Fic. 1. Block Diagram of Slow-Fast Coincidence Scintillation Spectrometer. 


Pulses from the anode of the photomultiplier are limited to an ampli- 
tude of 4 volt by a limiter circuit. The limiter output pulses are carried to 
the fast coincidence unit by a coaxial cable RG8U which is terminated at 
both ends by its characteristic impedance to eliminate any reflections. These 
pulses are shaped at the input to the coincidence circuit by 400 cm. long 
RG8U coaxial cable shorted at the other end. The fast coincidence circuit 
is a bridge,type circuit using crystal diodes. By inserting known delays in 
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either channel the prompt coincidence due to Co® gammas were counted and 
from the curve thus obtained, the resolving time of the coincidence unit was 
determined to be 1-5 x10- sec. 
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Pulses from the third dynode, after amplification, are fed to a single 
channel analyser for energy selection. The fast coincidence output and those 
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from the two single channel analysers are fed to a slow triple coincidence 
unit (resolving time 4 ys.) thus enabling coincidences between any two selected 
energies of gamma-rays to be counted. This indirect method of introducing 
energy selection to coincidence circuits is necessary, because the fast coinci- 
dence unit must be operated by the pulse derived from the first few photo- 
electrons produced in the photomultiplier, whereas the analyser is required 
to be triggered by a pulse proportional to the total number of photo-electrons 
released by the light scintillation from the phosphor. The analyser output 
pulse will therefore be considerably delayed with respect to the pulse operat- 
ing the fast coincidence unit and, moreover, this delay will vary with pulse 
amplitude. The fast coincidence circuit therefore cannot follow the analyser. 
The chance coincidences are also counted simultaneously by feeding the same 
three pulses to another triple coincidence unit (resolving time 4 us.) with one 
channel delayed by 10 us. 
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RESULTS 


Gamma-ray spectra in coincidence with 604 kev., 794 kev. and 1349 kev. 
respectively were taken. A very strong source was used. The following 
results are obtained as shown in Figs. 3, 4 and 5:— 
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Fic. 4. Spectrum in Coincidence with 604 kev. 

1. One channel was fixed to allow full photo-peak of 794 kev. as shown 
in the inset of Fig. 3. Both crystals were shielded with lead and were kept 
120° apart and were arranged to be 2” from the source. The gamma transi- 
tions in coincidence with the 794 kev. gamma-ray are found to be of energies 
604 + 10 and 555 + 15 kev. There is also strong evidence for the existence 
of a gamma-ray of energy 460 + 20 kev., which is not in coincidence with 
794 kev. (Fig. 3). 

2. One channel was fixed to allow a small fraction of the 604 kev. 
photo-peak as shown in the inset of Fig. 4. In this case the geometry was 
the same as above. The gamma-rays in coincidence with 604 kev. gamma-ray 
are found to be of energies 460 + 20 kev., 555 + 15 kev. and 794 + 15 kev. 

3. Finally, one channel was fixed to allow the full photo-peak of 
1349 kev. gamma-ray. Since this is weak, the crystal corresponding to this 
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channel was unshielded and kept at a distance of 1-5 cm. from the source. 
A lead plate, | cm. thick, was interposed between this crystal and the source 
in order to reduce the 604kev. and 794 kev. radiations. The 1349 key. 
gamma-ray is found to be in coincidence with 604 kev. gamma-ray but not 
with 794 kev. gamma-ray as seen from Fig. 5. 


DISCUSSION 


All the results obtained in this work, namely, the triple cascade 604 kev.- 
794 kev.—555 kev. and the double cascades 460 kev.—604 kev. and 1349 kev.- 
604 kev., are consistent with the decay-scheme of Cs!* shown in Fig. 6 which 
is the decay-scheme proposed by Forster and Wiggins." 


The main disagreement between the different decay-schemes proposed 
for Cs!*4 seems to be about the energy of the first excited state in Ba!™, 
Since the weak gamma 1349 kev. (1369 kev. of Forster!) is in cascade with 
the intense 604 kev. gamma-ray, 604 kev. as the energy of the first excited 
state seems to be more justified than 794 kev. If 794 kev. is assumed to be 
the energy of the first excited state, then the 1349 kev. gamma transition has 
to be to the ground state and the 604 kev. gamma-ray has to be emitted by 
the highest level 1972 kev. so as to be in cascade with the 1349 kev. transition. 
In that case, the 604 kev. gamma transition would be a weak one since the 
1972 kev. level is formed by a rather weak beta-ray group (80 kev., 28%). 
In fact 604 kev. gamma-ray has been found to be the most intense one. Hence 
it is preferable to assume 604 kev. as the energy of the first excited state. 


Author’s thanks are due to Dr. B. V. Thosar for his interest and helpful 
suggestions and to Dr. S. Jha for discussion. 
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FARADAY ROTATION IN CALCITE 
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Received August 23, 1956 


1. INTRODUCTION 


As the presence of birefringence in a solid interferes with the measurement 
of its magneto-optic rotation to a considerable degree, the determination of 
the Faraday rotation in a highly anisotropic crystal, even along the optic axis 
presents some experimental difficulty. In spite of its high birefringence, it 
was in calcite that Chauvin'»* verified the laws of propagation of light in a 
birefringent magneto-optic medium. His choice of this substance was ob- 
viously governed by the fact that perfect specimens of this crystal are easily 
available. Chauvin’s researches have demonstrated that the practical diffi- 
culties encountered in this field of experimentation are by no means insur- 
mountable. 


The magneto-optic rotation of calcite, apart from the problem it presents 
in its measurement, is interesting in many other respects. Previous researches 
from this laboratory on the magneto-optic properties of diamond,* the alums 
and sodium chlorate* have indicated that the presence of strong covalent bonds 
in the structure has a tendency to decrease the magneto-optic anomaly 
factor. Calcite with its CO, groups would, therefore, be an ideal substance 
to investigate. The present paper reports the values of the Faraday rotation 
in calcite parallel to the optic axis and also records the values of the Verdet 
constant of sodium chlorate. These measurements have proved quite fruit- 
ful in that it has been possible to establish, as had been foreseen by Darwin 
and Watson’ (1927) that the magneto-optic anomaly factor for various 
frequencies responsible for refractive dispersion may be different, 


2. THE VERDET CONSTANT IN CALCITE 


In a very good rhomb of calcite two faces were cut and polished per- 
pendicular to the optic axis. Great care was exercised (by a series of observa- 
tions under a polarisation microscope followed by grinding and polishing 
of the faces) to get the optic axis eaxctly normal to these planes. The soft- 
ness of the crystal, and the large variety of excellent abrasive and polishing 
papers that are available in the market to-day, considerably facilitated this 
process. The crystal was then fixed to a brass gadget (Chauvin?) which was 
capable of delicately rotating it through small angles about two perpendi- 
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cular axes. With this arrangement the crystal was set between the poles of 
the electromagnet so that light from a mercury arc passed exactly along the 
optic axis. The divergence of the light was reduced by first rendering it paral- 
lel and then allowing it to pass through four apertures, 3mm. in diameter, 
placed at intervals of 50cm. on either side of the calcite crystal. 


Measurement of the rotation was made for A 5461 with a sugar half shade 
(Rayleigh®), first placed at the polariser end and later at the analyser end. 
The magnitudes of the rotation with either arrangement did not differ by more 
than 4% showing that the birefringence and hence the missetting of the crystal 
was indeed small (Ramachandran and Ramaseshan’). The observations for 
A 5893, A 5780 and A5461 were made visually and those for A4358 and 
A 4047 were made by a spectroscopic technique developed by Landau,* de- 
tails of which have been given in an earlier publication. Table I gives the 
Verdet constant for calcite for different wavelengths. The specimen of cal- 
cite, 8-84 mm. thick, gave a rotation of 4-28° for a field of 16540 Oersteds. 


It is estimated that the absolute magnitude of the Verdet constants 
reported are accurate to 3% to 4%, while the errors in the relative dispersion 
values are less than 1%. 

TABLE I 
Magneto-Optic Data for Calcite 
V1 = Ye = 0-825, ys = 0-20. 





























Verdet Constant 
in min./cm./Oersted 
A y 
Exp. Calc. 
5893 | 0-0185 0:0184 0-488 
5790 | 00-0193 | 0-0192 0-488 | 
| 
| $461 00215 | 0-0216 0-481 
Be 
| 4358 00360 | 0:0355 .0°480 
| 4046 0-0418 0-0420 0-467 











3. THe MAGNETO-OpTiIC ANOMALY FACTOR 
From the measured values of the rotation and the known refractive dis- 


persion for the ordinary ray in calcite, the magneto-optic anomaly factor 
> was evaluated from the modified Becquerel formula 
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ae dn 
V=¥* gmc ** di @ 
These values have been entered in column 4 of Table I. One notices that 
with decrease in wavelength, y value gradually decreases from 0-488 to 0-467. 
As the anomaly factor is a measure of the deviation of the Zeeman splitting 
of the optical absorption frequency from the normal value, the variation of 7 
with wavelength can only signify that the anomaly factors associated with 
the various optical absorption frequencies must be different. 


If one assumes that the oscillator strengths (i.e., the transition probabi- 
lities) of the different dispersion frequencies are not affected by a magnetic 
field then if 


A,A? Api? 
e—l= y+ sm Gh atacucdarveisn (2) 





then it follows® that 


Be WrAA2A? —,yeAgAs2A? 7 
=n mat LO + | si 
Here y;, ya, ...., etc., are the anomaly factors for the absorption wavelengths 


Se 


In the case of calcite it is known that the refractive dispersion can be 
expressed by the formula 


_0-43257 #082930 =, 04337 2? 
A® — (0-0500)? ~ A®— (- 1000)? * A® — (0- 1535)? 
0-61855 A2 
+ 67? 


(Ramachandran’®) 


n,? —1= 





The magneto-optic dispersion is found to fit the formula 


a1 €, [ ,cronosezz as 0-006842 + 
~ m 2mc® L{A® — (0-0500)?}? © {A? — (0- 1000)?}* 
4, 0002044 12 
fr? — (0- 1535)" | 
from which y; = y2 = 0°825 and ys; = 0:20. Jn these formule V is in radians 


and X in microns. Column (3) of Table I gives the values of the Verdet con- 
stant calculated using these values of the magneto-optic anomaly factors. 
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The agreement between the experimental and calculated values is quite 
satisfactory. 


It is perhaps worthwhile at this stage to emphasise that the determination 
of y;, yo, etc., is actually by a “‘ fitting up” process. Therefore values of the 
anomaly factor of any accuracy or significance can be obtained only if one has 
(a) an accurate dispersion formula for refractive dispersion, using preferably 
experimentally observed absorption frequencies and also (b) very accurate 
values of the magneto-optic rotation over a wide range of wavelengths. In 
the case of calcite none of these conditions is very well satisfied. Even so, 
one notices that the anomaly factor for the absorption wavelength A 1535, 
which one could associate with the CO, group, is much less than that for the 
other wavelengths. 


4. FARADAY ROTATION IN NACLO, 

In view of this significant result it was considered worthwhile to redeter- 
mine the magneto-optic constants of NaClO, which also contains a group 
having strong covalent bonds. When the measurements were made it was 
found that the values reported by Voigt!® and later by Ramaseshan‘ contained 
a serious numerical error. The reported values were too small by a factor 
of two. Table II gives the Verdet constants for NaClOg,. 


TABLE II 
Magneto-Optic Data for NaClO, 


v1 = 0°74 ¥2 = 0°46. 





| | Verdet Constant 


| . ; , 
| in min./cm./Oersted 
| | 





Exp. | Calc. | 





5780 0)-0179 | 0-0179 0-668 | 

5461 0 0203 0-0203 | 0-669 | 

| 4358 0-0334 0-0334 0-663 | 

| 4047 0-0397 0:0397 | (0 660 | 

| 3650 0-0508 0-0507 0-653 | 

in this case ee Sede 
1-1825 A2 0-07992 A2 


~ 0:00864 A? 
(Chandrasekhar) 


nt—1l= 


A? — (0-0900)* ~ A? — (0-1850)? 
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and 


a 7 Tae lia =— (0: 7500) \2 T (2 — (0-1850)3}? 
from which y, = 0-74 and y. = 0°46. 





0-007088 A* . 0-001258 A* | 


The calculated values using these constants have been entered in column (3) 
while column (4) gives the values of the mean anomaly factor using formula (1). 


Here again one notices that the y values for an absorption wavelength 
associated with a covalent group (A 1850) is much lower than that for other 
wavelengths. It is, therefore, worthwhile to undertake the accurate measure- 
ment of the magneto-optic rotation, over a wide range of wavelengths for 
a large number of crystals. For such measurements a sensitive photo- 
electric polarimeter for visible and ultraviolet work has been developed by 
one of the authors (V.S.) and the results of the studies using it will be 
reported in due course. 

5. SUMMARY 


The Faraday rotation has been measured in calcite parallel to the optic 
axis for the wavelengths 5893, 5780, 15461, 14358 and ’4047A. A 
dispersion formula for the magnetic rotation has been fitted up using the same 
absorption frequencies that are used in refractive dispersion. It is shown 
that the magneto-optic anomaly factors for the various absorption frequencies 
are different and the value of the anomaly factor for the absorption wavelength 
A 1535, is very much lower than that for the other absorption wavelengths. 
(yYises = 0°20, yiooo = Ys00 = 0°825.) This paper also reports the values of 
the Verdet constant of NaClO, for a few wavelengths from A 5780 to A 3650 A. 
A numerical error in the earlier measurements has been noticed and rectified 
here. . In this case also a dispersion formula for magnetic rotation has been 
proposed and it is found that for this crystal y1959 = 0°46 while yo9, = 0°74. 
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(Communicated by Prof. R. S. Krishnan, F.a.sc.) 
1. INTRODUCTION 


ALTHOUGH crystalline quartz has been the subject of many investigations, 
very few accurate studies have been made on the magnetic and optical pro- 
perties of fused quartz. The dispersion of magnetic gyration in crystalline 
quartz has been measured by Borel (1903), Disch (1903), Lowry (1912) and 
Ingersoll (1917) over a wide range of wavelengths. But for fused quartz 
the Verdet constant has been determined for only three wavelengths A 5893, 
\ 5461 and A 4358 A (Stierlin, 1907; Lowry, 1912; Cotton, 1931; and Rama- 
seshan, 1946) and even here there seems to be considerable disagreement in 
the values reported. The refractive indices and the elastic constants of this 
substance also show variation from specimen to specimen. Till recently, most 
of the specimens of fused silica that were available, showed a curious strain 
structure probably caused by the structural birefringence arising from the 
previous thermal and mechanical history of the specimen (Raman, 1950). 
These flaws in the internal constitution of the specimens were obviously the 
cause for the discrepancies found in the measurements of the many physical 
properties of the substance. 

The Thermal Syndicate Ltd., London, has now been able to produce 
large specimens of very high optical quality fused quartz. A specimen of 
size 1” x 1” x 2” acquired by Prof. R. S. Krishnan from this concern for 
the study of light scattering and Raman effect was used in the present investi- 
gation. This specimen was found to be perfectly homogeneous and free from 
any inclusions, exhibiting a beautiful bluish track when a beam of light passed 
through it. Using this specimen, the absolute values of magnetic rota- 
tions for fused quartz have been determined accurately for different wave- 
lengths from A 5780 to 42848 A, and the results have been compared with 
those of crystalline quartz. The dispersion of magnetic gyration has also 
been correlated to the optical dispersion and the validity of the modified 
Becquerel formula has been tested. 


2. EXPERIMENTAL TECHNIQUE 


The magnet used in these experiments was of the Rutherford type fitted 
with special pole-pieces. Both visual and photoelectric techniques were em- 
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ployed for the determination of the magnetic rotations. A Lippich double field 
polariser and an analyser mounted on a divided circle that could read up 
to 0:01°, were used in the visual measurements. The complete experimental 
set-up and procedure for visual measurements have already been given in an 
earlier publication (Sivaramakrishnan, 1954). 


In the experimental set-up for the photoelectric method, the Lippich 
half-shade polariser was replaced by an ultraviolet nicol. The light from 
the analyser was focussed on to the slit of a monochromator capable of cover- 
ing the whole wavelength range in the visible and the ultraviolet. A R.C.A. 
1.P. 28 type photomultiplier tube was housed in a light-tight box attached 
to the exit slit of the monochromator. The power supply of 1000 V.D.C. 
for the photomultiplier tube was well smoothed and regulated. The com- 
plete wiring diagram of the photomultiplier with the power supply is given 
in Fig. 1. The output from the photomultiplier tube was either fed to 
a galvanometer directly or to a cathode-ray oscillograph through an A.C. 
amplifier. In the direct current measurement method, the dark current of 
the photomultiplier tube which is of the order of 10-* amps. was compensated 
by a counter electromotive force from a battery. This method was not how- 
ever extremely satisfactory when very great accuracy was aimed at, as the 
dark current was not constant but fluctuated erratically. 
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Fic. 1. Photomultiplier circuit diagram. 


Therefore for most: of the measurements, only the A.C. method of detec- 
tion was adopted, For this the incident light was modulated by running 
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the high pressure mercury discharge lamp source on the 50 cycles mains. The 
source was fairly stable in intensity as it had an intensity stabilizer unit. The 
signal from the photomultiplier was fed to a cathode-ray oscillograph through 
a narrow band amplifier. As in an A.C. amplifier the noise is proportional 
to the band width, the narrowing of the band width would considerably 
reduce the thermal noise from the phototubes, the hum from the power 
supply and the external sources which might arise due to imperfect shielding, 
etc. The amplifier used is similar to the one described by Miller, Long, 
Woodward and Thompson (1949) and Suryan (1953), but modified to suit 
the present investigation. The circuit diagram is given in Fig. 2. The 
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Fic. 2. Narrow band amplifier. 


amplifier was adjusted for maximum response at 100 cycles by the two filters 
F, and F, which were high pass and Twin-T network filters respectively. 
The amplification was controlled by means of three variable controls and 
extensive screening was done in the earlier stages to avoid the interference 
and spurious pick-ups. The sensitivity of detection could be increased further 
by using a phase detector with the narrow band amplifier. With such an 
arrangement, it would have been possible to detect rotations to a few seconds 
of an arc. However, various other factors limiting the accuracy of measure- 
ment, such as the stability of the magnetic field, the stability of the light 
intensity, and the maximum precision with which the analyser readings could 
be read, did not warrant the use of this more sensitive device. It was found 
possible, to measure the rotations accurately to -01° with the present set- 
up itself, 
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“The method of total extinction” and “the method of symmetrical 
angles” (Heller) were used for the measurement of rotations. With the 
same set-up it was possible to adopt both the methods. The experimental 
procedure is as follows. The substance under study was placed in between 
the pole-pieces of the magnet, and the two nicols were set for the total extinc- 
tion position. The analyser reading @ was noted. Later the magnetic field 
was put on and the analyser rotated so that there was minimum output from 
the photomultiplier tube as shown by the height of the peak of the 100 cycles 
sine wave in the cathode-ray oscillograph. Then the analyser reading 0, 
was noted and (@, — 4,) gave the approximate magnetic rotation for the 
substance for that particular field. Later the analyser was rotated and set 
at two positions A, and A, at which the output from the photomultiplier was 
equal, on either side of the approximate extinction position 6,. The mean 
of the two settings, i.e., (A, + A.)/2, gave the exact position of extinction. 
The rotations were measured with the magnetic field on and reversed. Usually 
the gain of the amplifier was altered by means of the attenuators and the 
rotations were determined at various amplifications. Also several readings 
were taken to avoid the errors due to fatigue effects. Measurements were 
made at number of wavelengths from 45780 to 42848 A and the rotations 
were found to be accurate to -01°. 


When solids with residual birefringence are placed in between the 
polariser and the analyser, it is not possible to set the two for toal extinction. 
This is overcome by compensating the signal from the photomultiplier tube 
(arising from the residual light intensity at the crossing position), by means 
of another signal from a second photomultiplier tube facing the same source 
of light. The two signals are mixed in a difference amplifier which precedes 
the main narrow band amplifier. 


The magnetic field was determined by measuring the rotation produced 
by conductivity water of known thickness (approximately the same thick- 
ness as the specimen). The Verdet constant for water at \5461A at 25° 
was taken to be -01547 minutes/cm./Oersted. (Rodger and Watson.) 


3. EXPERIMENTAL RESULTS 
The Verdet constant has been determined with light travelling perpendi- 
cular to all the three pairs of faces of the fused quartz block and it was found 
to be constant within 4%. The specimen was also rotated about an hori- 
zontal axis parallel to the direction of propagation of light and again the 
Verdet constant was practically constant with no systematic variation, thereby 
revealing the absence of any perceptible birefringence in the specimen, 
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(Ramachandran and Ramaseshan, 1952). The average value of the Verdet 
constant has been taken for the fitting up of a dispersion formula. 


In Table I are given the experimental values of the Verdet constant for 
different wavelengths. 


TABLE I 
Magneto-Optical Data for Fused Quartz 








a 




















| 
Verdet Constant | 
Wavelength Refractive | min./cm./Oersted | Mean 
A Index y value of 
Calc. Exp. | y 
ot —— 
| 
5780 1-4589 0-01538 0-01540 | 0-794 
5461 1°4601 0°01733 0:°01731 0-792 
4916 1 -4628 0°02173 0-02183 | 0-796 
4358 1 -4667 0-02828 0-02837 | 0-796 
4047 1 -4697 0-03337 0-03349 | 0°796 0-793 
3650 1-4744 0-04214 0-04224 | 0-796 | 
| 3342 1 -4800 0°0518 0-0516 | 0-790 | 
3132 1 -4845 0-06061 0-06059 0-793 | 
3022 1-4869 0-06614 0-06600 0-792 | 
2967 | -4887 0°06932 | 0-0690 0-790 | 
2848 1-4931 0-0768 | 0-0763 0-789 | 
| | 








It is worthwhile to compare the values of Verdet constant found in the 
present studies with those of Cotton (1931) and Ramaseshan (1946). The 
values of Verdet constant for the wavelengths A 5780, A 5461 and A 4358 of 
the above workers are given in Table II along with the values of the author. 


TABLE II 





| | Verdet constant min./cm./Oersted 
Wavelength A | 





Cotton Ramaseshan Author 











5780 0-01479 —-0-01491 0-01540 
5461 0°01671 001664 0-01731 
4358 002602 | 0-02640 0-02837 











It can be noticed that the values of Ramaseshan and Cotton differ from 
the author’s values by 4% at A 5461 and by 8% at A4358A, This may prob- 
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ably be due to the slight birefringence in the specimens studied by them. 
When birefringence is present, the observed rotations would be less than the 
actual ones, and the magnitude of the difference between the apparent and 
the actual rotations would increase with the decrease in wavelength (Siva- 
ramakrishnan and Ramaseshan, 1956). This is found to be true in the pre- 
sent case thus supporting the idea of the presence of residual birefringence 
in the specimens studied by previous authors. 


For evaluating the value of the anomaly factor we shall consider the 
dispersion formula proposed by Chandrasekhar (1951), 
0-5320 A? mie 0-4151 A? 0- 1570 A? 
12—(0-0600)2 * A®— (0- 1060)? * A®— (-1190) 
, Or40SB AF 1400 NT 
T NB (8-84)? " 2 — (20-74)? 


n—i1= 


which is found to fit the experimental data on refractive indices very satis- 
factorily. 

ey dn 
¥ Imc? “ dv 
V is the Verdet constant, y the weighted mean value of the anomaly factor 
of the different absorption frequencies, the value of y was evaluated at 
various wavelengths from the experimental values of Verdet constant and 
the total dispersion dn/dA obtained from the above formula for refraction. 
The evaluated values of » are given in Table I, Column 5. Itcan be noticed 
that the value of y for this substance in the wavelength region studied is prac- 
tically constant, the variation being less than 1%. This clearly indicates 
that all the absorption frequencies that are responsible to refractive disper- 
sion, are effective towards dispersion of magneto-optic rotation in the same 
proportion as in refraction. 


By using the modified Becquerel formula, V = where 


Therefore the following dispersion formula for magnetic gyration using 
all the absorption frequencies that are responsible towards refraction is 
proposed. 


© 1) g.993 {__0:0019150* | _—-004664.d 

~ Ime? n ~~ ~~ ((A® — (0-0600)?)2 © (A? — (0- 1060)?)? 
002223 a* i} 
" (A® — (0-1190)2}? 


Here V is in radians and A is in microns. 


The contribution to magnetic rotation by the infra-red frequencies is 
negligible and is not therefore taken into consideration. The value of the 











212 V. SIVARAMAKRISHNAN 


mean anomaly factor for fused quartz is 0-793. The Verdet constant for 
various wavelengths were calculated from the above formula, and are given in 
Table I. The calculated and the experimental values agree fairly well 
to 4%, which is the order of the experimental accuracy. 


4. MAGNETO-OPTIC ROTATION IN CRYSTALLINE AND FUSED QUARTZ 


As is well known, the formation of fused quartz from crystalline quartz 
when the latter is melted, is accompanied by a considerable diminution in 
the density from 2-651 to 2-203 and also diminution of the coefficient of 
thermal expansion which is 36 x 10-* for quartz and only 1-5 x 10-* for 
vitreous silica. We shall consider here the Faraday rotation and the optical 
refraction in these cases. The additive constant in the case of refractive 
index is refractivity and is given by 

n?—1 1 n? — | 

m+2a”" a 
depending on whether the Lorentz-Lorenz polarisation field is present 
inside the crystal or not. Similarly the additive constants for magnetic 
rotation is magnetic rotativity defined as 


ee a 
(n? + 2)? d d 
The significance of rotativity is that for a substance it remains constant even 


when the density changes. Table III gives the optical and the magneto- 
optical constants for crystalline and fused quartz. 











TABLE III 

‘ | | | i EI ear ae 
| A 5780 | Refractivity Rotativity 

| Refractive| Verdet | ee ey Oe) ieee ae are 

Substance Density index Constant | | | 
5780 | min./em./| 27-1 n*=—-j 1 nV | 9”  V{ 
_ Oersted ad n*+2 a d (2? +2)? d| 
a eee stead ERA: SAREE EAI ice | Sse 
Crystalline quartz 2-651 1-5447 *01714 | 0-5228 0-119) | *00999 -00467 
Fused quartz -+| 2-208 | 1-4590 | -01540 | 0-5111 | 0-1238 | 01017 | -00537 
| 


It can be noticed from the table that both the Verdet constant and the 
refractive index for fused silica are significantly lower than the corresponding 
constants for crystalline quartz. Also it can be seen that while there is a 
variation of only 2% in the values of nV/d for these substances, the variation 


i ° tn 
in Gzip g iS more than 10%. This is rather disturbing as _ these 
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results seem to throw a serious doubt on the very existence of Lorentz-Lorenz 
polarisation field in quartz. On the other hand this may also signify that 
fused quartz is fundamentally different in structure from crystalline quartz. 


We shall now consider the calculated values of anomaly factor for quartz 
from the values of Verdet constant reported by Borel. The optical dispersion 
dnidA was determined from the dispersion formula proposed by Radha- 
krishnan (1951), who makes use of the same absorption frequencies that are 
used in the refractive dispersion formula of fused quartz. Table IV gives 
Borel’s values of Verdet constant for crystalline quartz and the calculated 
value of y for various wavelengths. 


TABLE [IV 


Magneto-Optical Data for Crystalline Quartz 


: | er 
Wavelength | Refractive Verdet Constant : 
A Index min./cm./Oersted 
6438 «| s«1-5423 | 0-01368 | 0-773 
893 | 15443 | 0-01664 | 0-777 
5466 | 15462 | 0-01952 | 0-778 
5086 | 1 +5482 | 0-02257 | 0-773 
4799 =| ~—«1-5501 | 0-02574 | 0-778 
4678 15516 | 0-02750 0-785 
4358 | 1 -5538 | 0-03181 | 0-779 
3612 «| «11-5635 0-04617 | 0-740 
2573 =| = 1-5962 | 0-1079 | 0-747 
| | 








It can be seen from the table that the value of y drops down suddenly 
after 14358 A, for the two wavelengths for which measurements have been 
made. More extensive measurements in the ultraviolet region are needed 
before his data could be accepted as showing a real variation in y with wave- 
length. But if one considers only the values of Verdet constant from A 6438 A 
to 14358 A, the following dispersion formula fits the experimental data on 
magnetic gyration fairly well. 
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= 5,1 x 0-776 | 0002387 A? 0005819 A? 
~ Ime? n (A = (0-0600)2)2 ~ (A® — (0 1060)2)? 


, _0°002491 AP 
" (A® = (0-1190)?)25° 


The value of the anomaly factor 7 is 0-776. 


It is most interesting to note here that while the values of the Verdet 
constant of fused quartz is less than that of crystalline quartz, the anomaly 
factor y in the former is slightly greater than that in the latter. 


5. SUMMARY 


The complete experimental set-up of an accurate photoelectric polari- 
meter with an ultraviolet sensitive photomultiplier, a tuned amplifier followed 
by a cathode-ray oscillograph has been described in this paper. This has 
been used for the measurement of Faraday rotation in a perfect fused quartz 
specimen showing no residual birefringence. Measurements have been made 
for number of wavelengths A 5780 to A2848A. The following dispersion 
formula for magneto-optic rotation which fits the experimental data to $% 
has been proposed using the same absorption frequencies that are respon- 
sible for refractive dispersion, 





v=,£, ! x 0-793 [ O-OO1SIS.A® 0004664 A 
~ Ime? n (a2 — (0-0600)2)2 © (A? — (0-1060)2)2 
4 7_0'002223 A? 
QF =: ris07)"] 


where V is in radians and A in microns. 


The magneto-optic anomaly factor is 0-793 for all the absorption frequencies. 
The available data on crystalline quartz has also been analysed and the values 
of the Verdet constant, the magneto-optic anomaly factor and the magnetic 
rotativity of the two substances have been compared and found to be 
significantly different. 
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MAGNETO-OPTIC STUDIES IN THE NITRATES 
OF LEAD, BARIUM AND STRONTIUM 


By V. SIVARAMAKRISHNAN 
(Department of Physics, Indian Institute of Science, Bangalore-3) 
Received August 27, 1956 


(Communicated by Prof. R. S. Krishnan, F.A.sc.) 
1. INTRODUCTION 


MEASUREMENTS of the Faraday rotation in calcite and sodium chlorate have 
revealed that the magneto-optic anomaly factors for the different absorption 
frequencies of a substance may be different. It is also found that the anomaly 
factors for the frequencies which correspond to the anisotropic groups like 
CO;, ClO;, etc., are very low (Ramaseshan and Sivaramakrishnan, 1956). 
This idea has been supported by the studies of the author (1954) on ionic 
solutions where he found that the solvated nitrate ion had an extremely low 
value of 0-08 for the anomaly factor. The study of the magneto-optic pro- 
perties of solids with nitrate groups in them would therefore prove to be of 
great interest. The nitrates of lead, barium and strontium which were known 
to crystallise in the cubic system and could be quite easily obtained from 
aqueous solutions, were therefore chosen for the investigation reported in 
this paper. No measurements of the Faraday rotation in these crystals have 
so far been made. This is probably because of the fairly high residual birefrin- 
gence that these crystals exhibit. The method recently developed by 
Sivaramakrishnan and Ramaseshan (1956) has been employed here for eli- 
minating the effects of accidental birefringence in these crystals. The Verdet 
constant has been determined accurately for a number of wavelengths from 
A 5780 to A 3650 A using the photoelectric polarimetry technique. The values 
of the magneto-optic anomaly factor have been evaluated using the refractive 
dispersion values of Poinsot and Mathieu (1955). New dispersion formule 
for both optical refraction and magnetic gyration have been proposed here. 


2. MATERIALS AND METHODS 


The crystals of barium and strontium nitrates were grown in this labo- 
ratory by Dr. K. Vedam from warm-saturated solutions of the salts. Since 
the crystals were grown above 25° C., they were free from any water of hydra- 
tion. The specimens were 4 to 5 mm. in thickness and 2 sq. cm. in area and 
were found to be slightly birefringent. The crystal of lead nitrate used was 
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from Professor Krishnan’s collection and'\its size was 1-5 x | x lem. All 
these crystals exhibited a preferred axis of strain. 


The magneto-optic rotation in these crystals has been determined for 
a number of wavelengths from A 5780 to 43650 A. The measurements could 
not be extended below A 3600 A, since they were not appreciably transparent 
below that region. The photoelectric polarimetry method with A.C. com- 
pensation (Sivaramakrishnan, 1956) was adopted here for the measurements. 
The crystal was set such that the incident light vibration was parallel to the 
principal axis of strain. The apparent rotation , was then determined and 
later the polariser was rotated through 45° and ¥%,; was determined. From 
these values the true rotation p was evaluated using the formula 3 p = 2 #, + 
%,, (Sivaramakrishnan and Ramaseshan, 1956). Since the specimens of 
Ba(NO;), and Sr (NOs). were less than 0-5 cm. in thickness, the rotations 
were of the order of 1° to 2° only. Therefore the values of the Verdet con- 
stant for these two crystals are accurate to 1% only, while they are accurate 
to 4% in Pb(NOs).._ The effective magnetic field was estimated by measuring 
the rotation produced in conductivity water of the same thickness as the 
specimens. 


3. DISPERSION FORMULA FOR REFRACTION 


The data of Poinsot and Mathieu on the refractive indices of these crys- 
tals were used in the fitting up of the dispersion formule. Although the for- 
mulz proposed by these authors fit the experimental data on refractive 
indices well in the visible region of the spectrum, the fit in the ultraviolet 
region is very unsatisfactory. Therefore, new dispersion formule of the 
Drude type have been proposed using reasonable absorption frequencies. 


For the fitting up of a dispersion formula it is necessary to consider the 
absorption frequencies in the ultraviolet. But unfortunately the ultraviolet 
absorption frequencies have not been observed so far in these crystals. How- 
ever, it is known from the absorption spectra studies in aqueous solution of 
KNO,, that there is a strong absorption band near about A2100A. This 
has been ascribed to the nitrate ion (Maslakowez, 1928). Therefore, one of 
the absorption wavelengths for these crystals has been taken to be near about 
42100 A. Apart from this one, there may be many more absorption fre- 
quencies in these crystals. Since we do not know any of these frequencies 
experimentally, all these unknown frequencies are lumped up into one and 
dispersion formule having only two terms have been proposed here. More- 
over, since the data on refraction itself is accurate only to few units in the 
fourth decimal place, it is not worthwhile to attempt at a formula of many 
terms. 

Ab 
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The dispersion formule for refraction proposed now are as follows. 
Pb (NO), :— 


nt —1— 05740 146107, 09262 98 
~ 2® = (0-2250)? © A®— (0° 1048)? * AF — (7-143) 





Sr (NO,).:— 


pp OMNTE 1 24IT AE 
A* — (Q-2250)? ' A® — (0-1010)? 
Ba (NOs). — 


tp CRIT 1-768 A 
A* — (0-2250)? © A? — (0-0950)? 

It may be noted here that in these formule the complete contribution 
due to the infrared frequencies has not been taken into consideration. This 
is because of lack of data on refractive dispersion in the infrared and the 
little contribution made by these frequencies in the visible and the ultraviolet 
regions of the spectrum. Further we are interested only in the magneto- 
optic properties of these substances for which the infrared frequencies do not 
contribute much. However, it must be remembered that in order to satisfy 


the experimental data on dielectric constant, it is necessary to take them into 
consideration in the dispersion formule. 


The absorption wavelength 42250A which has been used for all the 
crystals may correspond to that of the nitrate ion. The second absorption 
wavelength varies with the substance. By using the above dispersion formule 
the refractive indices were calculated for various wavelengths and are given 
in Table I, along with the experimental ones and the difference between them. 


It can be noticed from the table that the experimental and the calculated 
values agree well to a few units in the fourth decimal place. This is quite 
satisfactory when compared with the formule proposed by Poinsot and 
Mathieu, where the difference between the calculated and the experimental 
values occurs in the second place of decimals for some wavelengths. 


4. MAGNETO-OpTiC ANOMALY FACTORS 


As pointed out by Ramaseshan and Sivaramakrishnan (1956), if we 
assume the refractive dispersion to be represented by a formula 


2 iA? 
then it follows that 
Vie | e y yiAjAj2r? 


~ n 2mc? ~ (A? — ),?)? 
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where yjs are the anomaly factors for the absorption wavelengths djs. If 
yiS are the same, then the value of the mean anomaly factor » in the modified 
Becquerel relation 
- @€ dn 

V=¥ omnct * Oh 
would be constant. From the dispersion formula for optical refraction pro- 
posed here, the values of dn/dA for various wavelengths were determined 
for these nitrates and the mean anomaly factors y were evaluated. It was 
found that the value of y decreased with the decrease in wavelength. The 
decrease in the value of y was found to b: 10% in Pb(NO,),. and more than 
20% in Ba (NOs). and Sr (NOs), crystals. This clearly indicates that the con- 
tributions to magnetic gyration by the two absorption wavelengths are not 
in the same proportion as in refraction. Therefore dispersion formule for 
magnetic gyration have been constructed with different values of anomaly 
factor for the two absorption wavelengths that are responsible for refractive 
dispersion. 


The formule are as follows: 








Pb (NO3)2:— 
y —¢ 1 ¢ 00103922 —_(0-01367 2? 
~ Zme* n ((A® — (0-2250)2)?  (A® — (0-1048)*)? 
Sr (NOs):— 
y — £1 § 0000658 A*_ 0-004979 2? 
~ Ime? n U(A® — (0-2250)2)? * (A® — O-1010)*) 
Ba (NO3)2:— 


Vn geay | fp OOM 5 O-OOSHIT | 
~ 2mc? n (A - (0-2250)?)? (A? — (0-0950)?)? 
where V is in radians and A in microns. 


The Verdet constants were calculated from these formule for various wave- 
lengths and are given in Table If along with the experimental values of 
Verdet constant and the evaluated values of mean anomaly factor y. The 
values of y, and y, for the three nitrates are as follows: 


Substance V1 Ye 
Lead nitrate es .. 0°36 0-85 
Strontium nitrate .. .. 0-061 0-40 


Barium nitrate aS .. 0°043 0:51 
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TABLE II 


Magneto-Optic Data for the Nitrates 


























Verdet Constant % | 
Substance Wavelength min./cm./Oersted y | 
In | 
cal. exp. | | 
| | 
lies ) | 
Pb (NO). . 5780 0-0492 0-0491 0-505 | 
5461 0-0563 0-0561 0-499 | 
4916 0-0732 0-0733 0-494 
4358 0-1010 - 0-1013 0-483 | 
4047 0-1253 0-1253 0-473 | 
3650 0-1743 0-1738 0-456 
Sr (NO,)s # 5780 | 0-0125 0-0124 0-222 
5461 0-0133 0-0133 0-216 
4358 0-0223 0-0225 | 0-202 
4047 0-0267 | 0:0264 | 0-189 
3650 0-0348 | 0-0348 | 0-177 
Ba (NO,), al 5780 0-0122 | 0-0123 {| 0-241 | 
5461 0-0138 | 0-0137 | 0-234 
4358 0-0228 | 0-0229 0-214 
4047 0-0270 | 0-0270 0-201 ‘| 
3650 0-0347 | — 0-0346 0-182 











| } | 





It can be easily seen from the table, that the experimental and the calculated 
values of Verdet constant agree fairly well within the limits of experimental 
error. It can also be seen that the value of the anomaly factor y,, correspond- 
ing to the absorption wavelength A 2250 A is very much lower than that for 
the other absorption wavelength in all the three crystals. In Ba(NO,). and 
Sr (NOs;)o, the values of (2259) are extremely low and are almost equal. In 
the studies on the ionic solutions by the author, it has been noticed that the 
value of the anomaly factor for the nitrate ion at infinite dilution is 0-08 
which strongly supports the view that the absoprtion wavelength at A 2250 
corresponds to that of the covalent nitrate ion. But in the case of Pb (NO,)., 
the value of the anomaly factor for the same absorption wavelength A 2250 A 
is higher and is equal to 0-36. This may possibly be due to the existence in 
the substances containing lead ions, of another absorption frequency near 
about 4 2000 A, whose y factor is very much higher. Studies of the author 
on the magneto-optic properties of lead glasses seem to lend support to this 
conjecture. 
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A critical study of the magnitude of the anomaly factor y and y. reveals 
some interesting features. It can be easily noticed that the value of > and y, 
in barium and strontium nitrates are lower than the corresponding values in 
the lead salt. This may probably be due to the tighter binding present in 
barium and strontium nitrates than in lead nitrate. Since a low value of the 
anomaly factor indicates a partial covalent binding, one could conclude that 
the binding in barium and strontium nitrates is more covalent than in lead 
nitrate. This conclusion is supported by other physical constants of these 
substances (Srinivasan, 1955). It is also supported from studies on Raman 
effect. The lattice lines appearing in the Raman spectrum of barium nitrate 
is found to be more intense than those in the lead salt (Couture & Mathieu, 
1947). 

5. SUMMARY 


The Verdet constant has been determined for the nitrates of lead, barium 
and strontium from 45780 to 13650A. Dispersion formule have been 
constructed for both optical refraction and magneto-optic rotation using the 
same absorption frequencies in both. In these nitrates it is found that the 
value of magneto-optic anomaly factor is different for the different absorption 
wavelengths of the substance and the value of y for the frequency corres- 
ponding to the nitrate ion is extremely low, and is equal to 0-06. The results 
indicate that the binding in barium and strontium nitrates may be more 
covalent than that in lead nitrate. 
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SYNTHESIS OF ANGOLENSIN 
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(Department of Chemistry, University of Delhi, Delhi-8) 
Received August 21, 1956 


ANGOLENSIN (Ia), an a-methyl-desoxy benzoin related to isoflavones, was 
first isolated by King, King and Warwick! from the wood of Pterocarpus 
angolensis and later by Gupta and Seshadri? from the wood and bark of 
Pterocarpus indicus. Its consitution (Ia) was established by degradative 
studies by King et al.’ but no synthetic support was provided. An obviously 
simple synthesis would be the C-methylation of the reactive methylene group 
of 2: 4-dihydroxy-4’-methoxy-desoxy benzoin. Badcock et al.* have shown 
that when heated with excess of methyl iodide and anhydrous potassium car- 
bonate in acetone solution 2-hydroxy-4-methoxy- and 2-hydroxy-4 : 6-di- 
methoxy desoxy benzoins suffer a-methylation besides O-methylation. But 
when there was a methoxyl group in the benzyl part no C-methylation took 
place. This reaction has now been examined in detail using various modifica- 
tions and has not been found to be suitable for the present synthesis. 


RO OH 
1 2 
i fn <S 2 
CH—. —OCH; CH;0— —CH=CH—CHy, 
+ \. \—/ 
I II 
a; R=H 
6; R=CHg 
3 SS “a. CHs - a. os Jens 
CH,0 2 CHC 0 CH,O . 3 ee 
Ill IV 


Consequently an alternative method involving Friedel and Craft’s reac- 
tion has been examined. p-Methoxy-hydratropic acid (IV) was needed for 
this purpose. This was originally prepared by Bougault* by the oxidation 
of anethole (II) first with iodine and mercuric oxide to p-methoxy-hydra- 
tropaldehyde (III) and then with silver oxide and alkali to p-methoxy hydra- 
tropic acid (IV). He reported that in many oxidations of the aldehyde to 
the acid the whole of the product got resinified. This has been our experience 
also. Hence the C-methylation of the ester and the nitrile of p-methoxy 
phenyl acetic acid was attempted with methyl iodide in the same way as with 
phenyl-acetonitrile and ethyl phenylacetate using sodamide® and also using 
sodium methoxide in methanol and anhydrous potssium carbonate in acetone 
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solution but the reactions were unsuccessful. This may be attributed to the 
adverse effect of the p-methoxyl group on the ionisation of the active methy- 
lene group. Finally the synthesis of p-methoxy-hydratropic acid (IV) has 
been conveniently carried out by the following route. 


Vv VI 


cuo-€ \>-co-cn, > emo _}=co-coon —> 


CHs 

; VA AN | . : foe 

CH,O— S-c—cooH > H,0— S—c 
ee . — _/ cooH 


OH 
VII VIII 

For the oxidation of p-methoxy-acetophenone (V) to p-methoxyphenyl 
glyoxalic acid (VI) with alkaline potassium permanganate a modification of 
the method of Bougault* has been employed giving a good yield. The second 
stage has been conveniently carried out by the Grignard reaction using methyl 
magnesium iodide in ether. p-Methoxy-atrolactinic acid (VII) was then 
dehydrated to p-methoxy-atropic acid (VIII) by following the procedure of 
Bougault’ and final reduction to p-methoxy-hydratropic acid (IV) was accom- 
plished using sodium amalgam and water. 


The feasibility of Friedel and Craft’s reaction using p-methoxy-hydra- 
tropic acid chloride was first tested using resorcinol dimethyl ether when 
angolensin-4-methyl ether (1 5) was obtained satisfactorily and it agreed in 
its properties with the earlier description of King et al. and with the sample 
prepared in this laboratory” from natural angolensin. Repetition of Friedel 
and Crafi’s reaction using O-dibenzyl resorcinol followed by debenzylation 
yielded angolensin (1 a) identical with the racemic form of angolensin isolated 
from the wood of Prerocarpus indicus.* 


—> IV 


EXPERIMENTAL 
p-Methoxy Phenyl Glyoxalic. Acid (VI) 

This acid was first obtained by Bougault® by alkaline potassium per- 
manganate oxidation of p-methoxy acetophenone in aqueous suspension at 
0°. Owing to sparing solubility of the ketone the yield is poor. The follow- 
ing modified procedure gives a much better yield. 


e 


A solution of potassium permanganate (32 g.) and potassium hydroxide 
(12 g.) in water (600 c.c.) was added dropwise to a solution of p-methoxy 
acetophenone (12 g.) in pyridine (100 c.c.) at 15° with stirring during a 
period of four hours. The stirring was continued for one hour more and 
the unreacted permangante decomposed with sodium bisulphite, Manganese 
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dioxide that had separated was filtered off and washed with two 100c.c. 
portions of water. The alkaline filtrate was acidified with concentrated 
hydrochloric acid (200 c.c.) and the solution left overnight in a refrigerator. 
Impure anisic acid (about | g.) that had separated out was filtered off, the 
filtrate saturated with sodium chloride, the solution extracted repeatedly with 
ether and the ether solution dried over magnesium sulphate. The residue 
after evaporation of ether crystallised from benzene-petroleum ether mixture 
as colourless thin and large rectangular plates, m.p. 87-88°. Yield, 6-7 g. 


p-Methoxy-Atrolactinic Acid (VII) 


This acid was obtained by Bougault? by permanganate oxidation of 
p-methoxy-hydratropic acid (IV). It has now been obtained from p-methoxy- 
phenyl glyoxalic acid (VI) by the Grignard reaction. 


Magnesium powder (1-6 g., 3 mols) was added to a solution of methyl 
iodide (4-2 c.c.; 3 mols) in dry ether (100 c.c.) containing a trace of iodine. 
When the metal had completely dissolved, the ether solution was cooled to 
0° and added dropwise to a solution of p-methoxy glyoxalic acid (4 g., 1 mol) 
in dry ether (100 c.c.) cooled in freezing mixture. After two hours standing 
in the freezing mixture a gummy solid was found to have separated out. The 
ether was decanted off and the solid treated with ice and dilute sulphuric acid 
when slowly it turned crystalline. It was extracted repeatedly with ether 
and the ether extract dried over magnesium sulphate. The residue after 
evaporation of ether crystallised from benzene as colourless prismatic needles 
and narrow rectangular plates, m.p. 129-30°. 


p-Methoxy-Atropic Acid (VIII) 


The following represents the best conditions for obtaining this acid by 
the method of Bougault.* A solution of p-methoxy-atrolactinic acid (5 g.) 
in glacial acetic acid (50 c.c.) was refluxed in an oil-bath at 130° for four hours. 
Acetic acid (about 40 c.c.) was distilled off under reduced pressure, the residue 
diluted with water (90 c.c.) and left at 0° for 24 hours when p-methoxy-atropic 
acid crystallised as colourless very thin plates, m.p. 118-20°. Yield, 2-5 g. 
A solution of the substance in glacial acetic acid decolorised bromine water. 
p-Methoxy-Hydratropic Acid (IV) 


Sodium amalgam (50 g., 2:5%) was added with constant shaking to a 
solution of p-methoxy-atropic acid (1 g.) in 2% aqueous sodium hydroxide 
(50 c.c.) and was left overnight. The brown alkaline solution was separated 
from mercury and acidified with concentrated hydrochloric acid when a semi- 
solid mass separated out which was repeatedly extracted with ether and the 
ether solution dried over magnesium sulphate, After the distillation of ether 
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the residue crystallised from benzene-petroleum ether mixture as colourless 
irregularly-shaped plates, m.p. 55-57°. Yield, 0-8 g. It did not decolorise 
bromine water in glacial acetic acid solution. 


Angolensin-4-methyl Ether (I b) 


The chloride of p-methoxy hydratropic acid (IV) was prepared as follows: 
To a solution of the acid (2 g.) in dry chloroform (25 c.c.) was added phos- 
phorus trichloride (0-6g., 1/3 mol.) dropwise with shaking and cooling. 
The reaction mixture was left at room temperature for two hours with occa- 
sional shaking, heated to 60° for two minutes and the chloroform solution 
separated from the phosphorous acid that had separated. On distilling off 
chloroform under reduced pressure the acid chloride was left behind as a 
very pale yellow viscous oil which was directly used for the Friedel and 
Craft’s reaction. 


A solution of anhydrous aluminium chloride (8 g.) was prepared in dry 
ether (80c.c.) with cooling in ice. To this was added resorcinol dimethyl 
ether (1-6 g.) and then an ether solution of the acid chloride prepared as 
above. The mixture was left at room temperature overnight. The next 
day as much of ether as possible was distilled off, the complex decomposed 
with ice and hydrochloric acid and extracted with ether. The ether solution 
was separated into aqueous sodium bicarbonate soluble fraction (A), sodium 
carbonate soluble fraction (B), sodium hydroxide soluble fraction (C) and 
the neutral fraction (D). Fraction (A) gave a small amount of unchanged 
acid; fraction (B) gave none on acidification and fraction (D) yielded a small 
amount of resorcinol dimethyl ether and hence all these were rejected. 
Fraction (C) on acidification gave a viscous oil which was extracted with 
ether and the ether solution dried over sodium sulphate. The residue left 
after removal of ether was distilled under reduced pressure when angolensin- 
4-methyl ether passed over at 230-33° at 2:3 mm. pressure. It gave a reddish 
brown colour with alcoholic ferric chloride and a green colour with concen- 
trated nitric acid. Yield, 1-2g. King et al.’ recorded 240-50° as bath 
temperature at 0-3 mm. pressure. 


The filter-paper chromatography of synthetic angolensin-4-methyl ether 
and the one prepared from natural angolensin by the procedure of Gupta 
and Seshadri? gave the same brown ring, Ry, 0°83 at 36° using 15% isopropyl 
alcohol as the solvent and diazotised benzidine reagent as the developer. The 
oxidation of synthetic angolensin-4-methyl ether with potassium permanga- 
nate in acetone solution was carried out by the procedure of Gupta and 
Seshadri? whereby p-methoxy-acetophenone, anisic acid and p-methoxy sali- 
cylic acid were isolated. 
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Angolensin (I a) 


Using the same procedure and resorcinol dibenzyl ether (3-2 g.) and the 
acid chloride from p-methoxy-hydratropic acid (2 g.) crude angolensin-4- 
benzyl ether was obtained. Yield, 0-5g. Small amounts of unchanged acid 
and resorcinol dibenzyl ether were recovered but the poor yield should be 
ascribed to the resinification that takes place during the reaction. 


The benzyl ether (0-5 g.) was dissolved in glacial acetic acid (20c.c.) 
and concentrated hydrochloric acid (20 c.c.) added. The mixture was heated 
on a boiling water-bath for one hour, poured over ice and extracted repeatedly 
with ether. The ether extract was washed with aqueous sodium bicarbonate 
to remove acetic acid and then extracted with 5% aqueous sodium carbonate. 
The carbonate solution was acidified, extracted with ether and the ether solu- 
tion dried over sodium sulphate. The residue after removal of ether crystal- 
lised from benzene-petroleum ether mixture as small prisms and plates, m.p. 
118-20° alone or admixed with an authentic sample of racemic angolensin 
isolated from Pterocarpus indicus.2 It gave a reddish brown colour with 
alcoholic ferric chloride, R;, 0-72 (brown ring) at 36° using the organic layer 
of methanol-chloroform-ligroin-water (1:2: 7:5) mixture as the solvent and 
diazotised benzidine reagent as the developer was obtained with both syn- 
thetic and natural samples and with a mixture of the two. 


SUMMARY 


By the Friedel and Craft’s reaction using resorcinol dimethyl ether and 
p-methoxy hydratropic acid chloride, angolensin monomethyl ether has been 
prepared. Employing resorcinol dibenzyl ether for the reaction and final 
debenzylation, angolensin (racemic) itself has been obtained. 
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INTRODUCTION 


FERNANDES! has reported the preparation of numerous compounds of various 
hydroxy organic acids with molybdic and tungstic acids and has mentioned 
three distinct complexes of molybdic oxide with potassium gallate in the ratio 
1, 2 and 0-66 moles of gallate per mole of MoO, to which the following 
formule have been assigned respectively :— 


MoO, C,H: (OH) CO;K; MoO, [O (OH).C.H:-CO,K],; MoO, —O——C,H, 
” oe 
VY 
MoO, 2 
(i) (ii) (iii) 


These were obtained by treating solutions of potassium gallate with 
calculated quantities of the molybdenum oxide and concentrating the solu- 
tions when the compounds crystallised out. However in a subsequent com- 
munication? on the complexes of other hydroxyacids with MoO,, WO, and 
UO, it has also been stated that in aqueous solutions compounds of the 
first type, viz., 1: 1 complexes, are normally obtained but a very large excess 
of the phenol is required to give complexes of the second type. 


Pozzi-Escot* has utilised the yellow to orange colour given by molyb- 
dates in presence of polyphenols, such as gallic and pyrogallic acids, for the 
detection of molybdenum. 


The work presented in this paper was undertaken firstly to find out by 
spectrophotometric method, under controlled conditions of pH, etc., the 
evidence for the formation of the various compounds reported and secondly 
to find out the feasibility of utilising the colour formation for the quanti- 
tative estimation of molybdenum, 


228 

















Studies of Gallic Acid Complexes with Metals—II 


EXPERIMENTAL 
Apparatus and Reagents 


Instruments.—pH measurements were made on a Marconi pH meter, 
type TF 511 D. Spectrophotometric measurements were done on a Beckman 
Model DU Quartz Spectrophotometer using | cm. quartz absorption cells. 


Gallic acid —E. Merck, G. R. recrystallised and dried at 120°C. Fresh 
solutions were prepared when required. 


Sodium molybdate-—B.D.H. Analar molybdic acid dissolved in a little 
more than the stoichiometric quantity of sodium hydroxide solution (from 
E. Merck, G.R. grade) and standardised by gravimetric estimation as lead 
molybdate. 


Sodium acetate.—E. Merck, G.R. grade. 
Sulphuric acid —E. Merck, G.R. grade. 


ABSORPTION SPECTRA 


The absorption spectra of the two reactants, viz., sodium molybdate 
in 2-8 x 10° M concentration, gallic acid in 1-5 x 10-*M concentration and 
of the complex formed by mixing these, all at a pH of 4-2, are presented in 
Fig. 1. The optical density values for the complex have been obtained after 
allowing for the absorption due to the two reactants. The curve shows a 
single absorption peak at 305 my. 


Absorption measurements were also carried out with more concen- 
trated solutions in the region of 300-600 my using ratios of molybdic acid: 
gallic acid as 1:1, 1:4 and 1:10°5. No characteristic absorption was 
observed in any case as seen from Fig. 2 (plotted up to 500 my only). 


EFFECT OF PH 


Preliminary experiments had shown that the formation of the molyb- 
denum-gallic acid colour was extremely sensitive to hydrogen-ion concentra- 
tion and showed maximum absorption in the pH range of 4-4-5 beyond 
which a sharp decrease in the absorption was noticed. The results obtained 
at 320 mp are shown in Fig. 3. As the absorption of gallic acid is affected 
by the change in pH,*: ® this effect was obviated by taking gallic acid solutions 
of corresponding pH as blanks. 


The same effect was observed in the absorption spectra in the region 
of 250-375 mu by taking molybdic acid and gallic acid in the ratio of 1:4 
and changing the pH (Fig. 4). 
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Fic. 1. Absorption Spectra:—Curve I: 2-8x10-°M. Molybdate. 
» Il: 1-5x10-*M. Gallic Acid. 


5, II: Complex, 2°8x10-°M. in Molybdate and 
1-5 x10-* M. in Gallic Acid (all at pH 4-2) 


In all further experiments therefore the pH was adjusted to 4-2. 
COMPOSITION OF THE COMPLEX 


The ratio of Mo to gallic acid in the complex or complexes was investi- 
gated by three different methods, viz., (a) method of continuous variations,°*’ 
(6) the molar ratio method® and (c) the modified slope ratio method.®:” 
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Since molybdenum is stated to form more than one complex with gallic 
acid it was necessary to find out suitable wavelengths at which measurements 
for the various complexes could be made. 


Absorption spectra of the solutions were first obtained where molar 
ratios of Mo to gallic acid were varied. These absorption curves above the 
wavelength of 300 my are shown in Fig. 2. Owing to high absorption of 
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gallic acid in the ultraviolet region it was not possible to make measurements 
in the region lower than 300 my for solutions containing high ratio of gallic 
acid to molybdenum. Similarly curves ll and IV of Fig. 4 show the absorp- 
tion spectra of solutions containing molybdenum: gallic acid in the ratios of 
1:4 and 1:1 respectively all the other conditions being constant. 
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The spectral data does not in the least indicate the formation of more 
than one complex in the solutions or alternately if more complexes do exist 
they have similar absorption characteristics. 


(a) Method of Continuous Variations 


In the absence of any characteristic wavelengths of absorption for differ- 
ent complexes, if at all they are present, the study of the system by the method 
of continuous variation was made at three wavelengths, viz., 300 and 315 my 
which are near the maximum absorption wavelength of the complex and at 
360 mp where the reactants themselves will have a very low absorption. 


A series of solutions was prepared in which the ratios of molybdenum 
to gallic acid was varied but the total molarity of the two was constant, the 
PH also being kept constant at 4-2, and their optical densities were measured. 
f Fig. 5 shows the results obtained at wavelengths of 300, 315 and 360 mp 
A6 
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Fic. 5. Continuous Variations. 


where the difference (D) between the observed optical density of each solu- 
tion and that calculated on the assumption that there is no reaction between 
the two reagents have been plotted against the ratio [Mo]/[Mo] + [C,H,0,]. 
For convenience of representation the values of D for the wavelength of 
300 my have been multiplied by a constant factor in Fig. 5. 


In all the three curves the maximum absorption occurs when the relative 
proportion of gallic acid to molybdenum is 1:1 forming a complex which 
from the nature of the graphs indicates that it is highly dissociating in solu- 
tion. 


(b) Molar Ratio Method 


Fig. 6 shows the effect of increasing the molar concentration of one of 
the reagents at a fixed concentration of the other, the measurements being 
made at 380 mp. Owing to the high dissociation of the complex the curves 
do not show any sharp break. Attempts were made to get a sharp break in 
the curve by increasing the ionic concentration of the solutions by addition 








Studies of Gallic Acid Complexes with Metals—II 235 
A slight increase in the optical densities was observed 


of potassium chloride. 
but the increase was so slight as not to be of any consequence. 
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the point of intersection corresponds to a molar ratio of molybdenum: gallic 
acid as 1:1. 


Further it may be noted that curve I (Fig. 6) in which the concentration 
of molybdenum is varied with respect to a constant concentration of galiic 
acid levels off to a reasonably constant optical density, whereas in curve II 
the optical density was still found to be rising even at a molybdenum: gallic 
acid ratio of 1:6. In fact increasing the ratio of molybdenum: gallic acid 
up to 1:50 still showed increasing values in optical density. This may 
possibly be due to the formation of only 1 : 1 complex in the former case while 
in the latter case the second complex of molybdenum (Mo: gallic acid:: 1:2 
as previously reported) is probably formed as more and more gallic acid is 
made available. No direct evidence for the formation of the molybdenum: 
gallic acid:: 1: 2 complex was however obtained. 


(c) Modified Slope Ratio Method 


Raghava Rao and co-workers?, #° have used a modified slope ratio method 
in which they plot the optical densities versus the concentration of the 
variable components from either end of the X-axis, thus giving two plots 
which intersect at a point corresponding to the composition of the complex, 
the assumption being that only one complex is formed. 


In the discussion under “ Molar ratio method ” it has been shown that 
a large excess of gallic acid increases the optical density of the solutions con- 
tinuously. Restricting therefore to molar ratio of one reactant up to 2 with 
respect to 1 mole of the other a plot is obtained as given in Fig. 7, which once 
again gives the ratio of molybdenum: gallic acid in the complex as 1: 1. 


DISSOCIATION CONSTANT 


The dissociation constant of the complex was calculated by the method 
described by Anderson and co-workers." 12 Assuming that the extreme end 
of curve I of Fig. 6 represents complete colour saturation for the quantity 
of gallic acid used (0-0005 gramme moles in one litre), the molar extinction 
coefficient of the complex was calculated. From the absorbancy values and 
the known concentrations of the two reactants, the equilibrium concentra- 
tions of the complex and the two reactants were calculated. From these 
the dissociation constant (K) was calculated. The average value for 

_— [Mo] [R] 


where R represents gallic acid. 
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Fic. 7. Modified Slope Ratio. 


ESTIMATION OF MOLYBDENUM 


Contrary to the observations of Pozzi-Escot® it was found that by main- 
taining pH constant at 4-2 Beer’s law was valid for the ranges studied, viz., 
0-01 to 0-10 mg. and 0-10 to 1 mg. of molybdenum per 100 c.c. of solution 
in the presence of 15mg. and 150 mg. respectively of gallic acid. 
Measurements were preferably taken at wavelengths of 330 mp and above as 
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the absorbancy due to gallic acid is low in this region. The colours are 
developed immediately and were found to be stable up to 24 hours. 


The graphs as given in Fig. 8 for the two ranges of molybdenum con- 
centration studied show different slopes. It may be remembered however 
that the excess of gallic acid used in the two sets is very much different. 
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INTERFERENCES 


A number of metal ions such as Fe, Ti, V, W, etc., interfered in the esti- 
mation. This interference could not be easily removed. Addition of com- 
plexing agents like tartaric and oxalic acid or E.D.T.A. either diminished the 
molybdenum colour or destroyed it completely. It is therefore necessary to 
isolate the molyebdnum before developing the colour. 


SENSITIVITY 


Rice and Yerkes!* have reported the sensitivities of various reagents 
used for molybdenum estimation. A comparison of these with the present 
method is given below :— 


Potassium ferrocyanide = ‘ia 7” 0-1 
Sodium thiosulphate . . ius dan ‘i 0-04 
Phenyl hydrazine si oe ia ‘. 0-006 
Potassium xanthate .. ree ea “A 0-004 
Butyl xanthate <n ee = ie 0-004 
Amyl xanthate - i re - 0-004 
Pentasol xanthate .. na ‘ 0-004 


Potassium ethyl xanthate (colour extracted with 0-0004 
ether) 


Gallic acid .. - i ‘i a 0-0001 (O.D. = 0-01) 
Ammonium thiocyanate-stannous chloride (colour 0-00002 
extracted with ether) 


SUMMARY 


Spectrophotometric study of the molybdenum-gallic acid system showed 
the formation of a yellow complex of the empirical formula MoR which shows 
an absorption maximum at 305mp. The complex dissociates appreciably 
in solution and gave the value for K, the dissociation constant as 0:6 x 10~ 
+0-1 x 10. Beer’s law was found to hold good for concentrations up to 
1 mg. Mo per 100c.c. and quantitative estimations could be made after iso- 
lation of Mo from interfering elements. 


The authors are thankful to Mr. M. S. Das, Mr. Ch. Venkateswarlu 
and Mr. G. K. Belekar for valuable suggestions. 
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INTRODUCTION 


THE reaction between gallic acid and ferrous sulphate in tartrate medium 
has been employed for the colorimetric estimation of gallic acid.!: 2» > Glass- 
tone* has reported the use of osmium tetroxide as a colorimetric reagent. 
Mattil and Filer? have developed a spectrophotometric method based on the 
ultraviolet absorption of gallic acid itself at 270-71 my in 0-01 N hydrochloric 
acid. 


In connection with some work on the metal-gallic acid complexes in this 
laboratory (yet to be communicated), a rapid colorimetric method for esti- 
mating gallic acid was required. For this purpose, none of the above- 
mentioned methods were found to be very well suited chiefly because of the 
very low concentration of gallic acid (1-6 mg. gallic acid/100 ml.—maximum) 
estimated. 


It was felt that the titanium-gallic acid system® could be utilised for the 
estimation of gallic acid, provided the hydrolysis of the excess titanium salt, 
required for colour saturation, could be prevented by the addition of another 
complexing agent, which would not otherwise interfere with the formation 
of the coloured complex. To this end, the effect of tartrate, oxalate and citrate 
ions on this system has been studied and the results are described in this 
paper. 

Apparatus and Reagents 


Same as in Part I (loc. cit.) 


Ammonium oxalate ) 
Ammonium tartrate + B.D.H., AnalaR 


Ammonium citrate J 


EXPERIMENTAL 


Preliminary experiments are carried out using ammonium tartrate as 
the complexing agent for the excess titanium. The relative proportion of 
ammonium tartrate to titania is arbitrarily fixed as 10: 1 (w/w) in the begin- 
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ning. The pH of the final solutions is always about 5-0. Under these condi- 
tions, the following observations are made: 


(i) The coloured complex shows maximum absorption at 370 my in 
the region 300-500 my. 


(ii) The colorimetric titration at 370 mp with 2-0 mg. of gallic acid in 
a total volume of 100 ml. reveals that, though the colour saturation is more 
or less complete when 2-5 mg. of titania and 25 mg. of ammonium tartrate 
are added, it is not quantitative even after the addition of 10 mg. of titania 
and 100 mg. of the tartrate. 


(iii) The calibration curves, studied between 370 and 450 my are linear 
showing that Beer’s Law is obeyed but they show a varying positive Y inter- 
cept and, 


(iv) the extinction values for the same amount of gallic acid are varying 
with variation of the quantity of titania used as well as with the relative pro- 
portion of the tartrate to titania. 


In view of the above observations, the role of tartrate in the system is 
further investigated as follows :— 


1-0 mg. of titania and 0-25 mg. gallic acid (mole ratio: approximately 
8: 1) are taken in different flasks. Varying amounts of the tartrate are added, 
pH adjusted to about 5-0 with 2 N ammonium acetate, and the volumes 
made up to 50 ml. Extinctions are measured at 370 my and plotted against 
the concentration of the tartrate (Fig. 1 a). It will be seen that the extinction 
first goes down, then rises to a more or less steady value and drops again 
slowly with further addition of tartrate. It is evident from this, that colour 
development is maximum when the tartrate concentration is 20 to 28 mg. 
(approximately corresponds to 9-12 moles of the tartrate for a mole of titania). 
The same experiment is repeated with different amounts of titania and the 
same amount of gallic acid, and in every case, the maximum colour formation 
is obtained when the ratio of the tartrate to titania is the same as found before. 


Similar experiments using ammonium oxalate and citrate (Figs. 1 b and 
1 c) show that these two complexing agents do prevent the hydrolysis of the 
titania but interfere seriously with the formation of the coloured complex. 


In all further experiments, a ratio of 1:25 (w/w) between titania and 
ammonium tartrate is maintained and this mixture is hereafter referred to 
as ‘titanium tartrate,’ the concentration of the mixture being expressed in 
terms of the titania content. 
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Absorption Spectrum of the Complex 


0-5 mg. of gallic acid is complexed with excess ‘ titanium tartrate’ in a 
total volume of 50 ml. and the absorption spectrum has been studied against 
reagent blank between 300 and 500 my. It shows maximum absorption at 
370 mp (Fig. 2). 


Effect of pH 


The effect of hydrogen-ion concentration is studied in the presence of 
tartrate. Maximum colour development is attained at pH 4-0-7-0 (same as 
in the case of TiGa,5). 
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Fic. 2. Absorption spectrum of the complex with 0-25 mg. gallic acid in 50 ml. 





Effect of Reagent Concentration 

The colorimetric titration is carried out with 1-0 mg. gallic acid in 50 ml. 
following the same procedure given in Part I (Joc. cit.). It is evident from 
Fig. 3 that colour development reaches the saturation limit this time, whereas 
preliminary experiments with a different amount of the tartrate have shown 
a continuous increase with increase of the reagent. 


Estimation of Gallic Acid 


Varying amounts of gallic acid are taken, ‘ titanium tartrate’ is added 
in excess (2:5 to 5-0 mg. of TiO, per 1-0 mg. of gallic acid), pH is adjusted 
as before with 2 N ammonium acetate to about 5-0 and the volumes are made 
up to 50 ml. The optical densities are measured at 370 my or at higher wave- 
length depending on the concentration of gallic acid taken, and plotted 
against the concentration of gallic acid. The slopes for a concentration range 
of 0-1 to 2-0 mg. of gallic acid in 50 ml. at 370, 400 and 450 mp are shown 
in Fig. 4. Now the.calibration}curves are not only linear but also pass through 
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Fic. 4. Calibration curves. 
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the origin. Using the above procedure, the maximum amount of gallic acid 
estimated is 10 mg. in a total volume of 50 ml. The method, thus, is accurate 
and very rapid overea wide range of gallic acid concentration. 


For the calibration purposes, only ‘ M and B’” titania is used. 


A study of the entire titanium-gallic acid system in the presence of these 
complexing agents is under further investigation. 


SUMMARY 


A comparative study of the influence of tartrate, oxalate and citrate 
ions on the titanium-gallic acid system, when titania is in excess, shows that 
the oxalate and citrate ions interfere seriously in the formation of the coloured 
complex. In the case of tartrate, maximum colour development is attained 
when the amount of ammonium tartrate is 20 to 28 times (w/w) to that of 
titania. Based on this, a spectrophotometric method has been described 
for the estimation of gallic acid over a wide concentration range. 
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